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Abstract: 

High energy scattering processes involving jets generically involve matrix elements of light- 
like Wilson lines, known as soft functions. These describe the structure of soft contributions 
to observables and encode color and kinematic correlations between jets. We compute the 
dijet soft function to O(a^) as a function of the two jet invariant masses, focusing on terms 
that have a non-separable dependence on these masses and are not determined by the 
renormalization group evolution of the soft function. Our results include non-global single 
and double logarithms, and analytic results for the full set of non-logarithmic contributions 
as well. Using a recent result for the thrust constant, we present the complete 0{a1) soft 
function for dijet production in both position and momentum space. 
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1. Introduction 



Jets play a ubiquitous role at all high-energy colliders, both as signals for new particles 
which interact strongly, and as backgrounds for such signals through Standard Model 
processes involving the strong interactions. Thus it is crucial to achieve reliable and precise 
theoretical predictions for many types of jet cross sections. 

Two obstacles to achieving precision predictions for jet observables are poor behavior 
of perturbative expansions and uncontrolled nonperturbative corrections. Two key tools 
for overcoming these challenges are factorization and resummation of large logarithms Q . 
For a jet cross section which factorizes, the separation of perturbative and nonperturbative 
corrections increases the predictive power thanks to universality of the nonperturbative 
contributions, and the perturbative contributions can be organized according to the differ- 
ent dynamical scales (e.g. hard, jet, soft) contributing to the cross section. The presence 
of these disparate scales allows factorization but is also the culprit in producing large loga- 
rithms in perturbative results for jet cross sections. Techniques to resum these logarithms 
to all orders in perturbation theory are often critical to obtaining well-behaved predictions. 

The most powerful methods to resum such logarithms rely on the renormalization group 
evolution of the different factors that appear in a factorized jet cross section, which are 
separated at an arbitrary scale fi. Such methods have been developed directly within the 
language of perturbative QCD itself ||, ||. An alternative approach has been formulated 
in the powerful language of the soft-collinear effective field theory (SCET) []|, |6|, [7|, ||, ||] . For 
e + e~ annihilation cross sections in which the jet-like structure of an entire event is probed 
with a single ("global") variable, such as thrust [ 10 1 , heavy jet mass [11], angularities 



12], etc., these methods succeed in resumming all logarithms of the "event shape" that 



become large in the two-jet limit to a well defined order in resummed perturbation theory 
H [3|, [l3|, 0, [l5|, |l6], [17|, [L8). For such a global two-jet event shape, e, the cross section 
takes the form da/de = da/de\ di . ct [l + 0(e)], where the leading order piece factorizes 
schematically as |l2|, jljj, 

da 

de 



= aoH(Q;v)[j n (Qe 1 / j ;fx)®J n (Qe 1 / j ;v)®S(Qe-,v)] , (1.1) 

dijet 



in terms of a hard function H, jet functions J n> n, and soft function S. The (g) denotes 
convolutions in the e-dependent arguments of J n ^ and S. For event shapes like thrust and 
jet mass, the exponent j = 2, but for angularities, j takes a range of values greater than 1. 
<7o is the total Born cross section. Each function depends on logs only of a single ratio of 
scales, n/Q for H, ^/(Qe 1 ^) for J, and ///(Qe) for S. Solving the renormalization group 
evolution equations in \i for each of these functions produces a form for da/de in which 
all large logarithms of e are resummed to a given order in resummed perturbation theory. 
The same technology can also be used to factorize and resum cumulative distributions or 
"cumulant" observables, 

As a step towards more exclusive probes of jets than Eq. (|El]), we can consider cross 
sections differential in more than one measure of the "jettiness" of a final state. We 
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will focus here on the particular example of dijet cross sections in the context of e + e~ 
annihilation at center-of-mass energy Q, in particular, the dijet invariant mass distribution 
to hadronic final states X defined by 



x 5^m\- (j^Pi^j j<^™|- 




(1.3) 



where and L M are sums of the vector and axial currents in QCD and QED, respectively. 
L, R are the two hemispheres defined with respect to the thrust axis of the final state X, 
and pi is the four-momentum of the ith particle in L, R. Using the formalism of SCET 
H> H> @> & [J> ^ nas been shown that this cross section factorizes into the form 13, 15], 



a H(Q,fi) I M 1 dl 2 J n {m\-Qt 1 ,ii)J n {ml-Qi 2 ,ii)S{h,l2^) + --- , (1.4) 



dm 2 dm 2 



to all orders in a s . The soft function S(£\,£ 2 , n) is given by a matrix element of Wilson 
lines and is defined below in Eq. (|2.1|) . The ellipses denote that the result is at leading 
order in the power expansion in m 2 2 /Q 2 <C 1. The corresponding factorization theorem 
for the double cumulant S(m^ 2 , mf) = / dm\dm\ 9(m^ — rn^Qim^ — m^) d 2 a j (dm^dm^) 
is 

E(mf ,mf) = a H(Q, f i)Q 2 Jd£ c 1 d£\ J n (mf - Ql\ , /i) J n (mf - Qi\ , ^)S C (l\,t 2: ix) + • • • , 

(1.5) 

where S r is the double cumulant soft function 



1 



S c {e c i,£h»)= I dh I d£ 2 S(£ 1 ,£ 2 ,»). (1-6) 



Eqs. (|1.4j ) and ( |1.5| ) exhibit a richer structure than Eq. (IT); containing more informa- 
tion about the two-jet-like final state, and presenting additional challenges to resumming 
all potentially large logarithms. Each of the hard and jet functions depend only on ratios 
of [i to a single scale, Q or rrii, but in contrast to Eq. ( |1.1|) the soft function now depends 
on the ratios fi/£i, {J-/£ 2 , and £\/£ 2 . When m\ ~ m 2 the factorization theorem Eqs. ( |1.4| ) 
and ( |1.5|) allow resummation of logs of mi /Q and m 2 /Q to arbitrarily high accuracy. Since 
it leaves logs of m\jm 2 in fixed-order perturbation theory without resummation, it does 
not handle log resummation for m\ 3> m 2 . In the latter situation hierarchies appear in the 
ratios appearing in the hemisphere soft function. These latter logs are examples of what 
have been dubbed "non-global logarithms" (NGLs) |2l[ . 

At present, all of the ingredients in Eq. (|l]4|) are known analytically at 0{a 2 s ) except 
for the soft function S(£i,£ 2 , fi). What is known so far about S at this order are all the logs 
of iij£\ and \ij£ 2 thanks to knowledge of the anomalous dimensions to two- loop order (the 
cusp anomalous dimension is known to three- loop order). These terms are constrained by 
renormalization group invariance of the physical cross section Eq. ( |1.4j ) , requiring the sum 



-3- 



of the anomalous dimensions of H, J n ,Jn, and S to add to zero. What is not known are the 
functions of h/h that can arise in S at 0(a 2 ). Conjectures have been made about what 



types of functions can arise [g2|, but these have never been validated nor their coefficients 
calculated analytically. 

Knowledge of the soft function S(l\,l2,p) is applicable not only in the dijet invari- 
ant mass distribution Eq. ( |1.4| ) itself, but also to a wide class of event shapes in e + e~ 
annihilation, including thrust, heavy jet mass, and the "asymmetric" thrust and jet mass 
variables [22|. Knowledge of S to 0{a 2 s ) is required to achieve NNLL (and higher) accuracy 



in resummed predictions for distributions in these observables, and is important input for 



recent extractions of the strong coupling a s at N 3 LL from thrust [16, n8f and the heavy 



jet mass f23||. These analyses currently depend on numerical extraction of 0(a 2 ) constants 



from the Monte Carlo generator EVENT2 [24, Pq| . The same S also enters event shape dis 



tributions for massive jets fl3| , ffq , E6j, applicable for example for extracting the top quark 



mass from jet mass distributions. Although Eq. (1.4) is formulated for e + e~ collisions, the 
hemisphere dijet soft function is actually closely related to an incoming dijet soft function 
that appears in event shapes for hadron-hadron collisions, such as the "beam thrust" or 
"O-jettiness" cross sections introduced in p7j , E8|. In that case, the two masses are those 
of the measured radiation in hemispheres determined by the beam directions. 

Observables like the dijet invariant mass distribution can probe the jet-like structure 
of an event in a "non-global" way, meaning that they are sensitive to soft radiation at 
different scales in sharply divided regions of phase space. The remaining sensitivity to 
these scales in the soft functions produces NGLs. Other examples of such non-global 



observables are exclusive jet cross sections [29] and jet shape distributions [[3(J, 31 1 using 



particular jet algorithms, as well as distributions from jet substructure algorithms p2| . 
Accounting for and resumming these logarithms will be important for achieving precision 
jet phenomenology in this era where we probe jets with ever more exclusive measures. 

NGLs thus far have not been resummed using the renormalization group-based tech- 
niques mentioned above. Their presence was first pointed out by Dasgupta and Salam in 



|2C, 21]. They arise in observables probing soft radiation only in a part of phase space (type 
1) or in sharply divided parts of phase space with different scales (type 2). An example 
of the first type was given in |20||, which considered the cumulative sing/e-hemisphere in- 



variant mass distribution E(/9r), related to our dijet invariant mass distribution Eq. (1.4) 
by 

da 



poo p 

= / dm i 

Jo Jo 



2 dm 2 dm 2 



Dasgupta and Salam noted (in our language 1 ) that a prediction for X(/°b) based on inserting 
the dijet factorization theorem Eq. (|1.4| ) in the relation Eq. (|1.7| ) would not resum all 
logs of pr from renormalization group evolution. Resummed results therefore need to be 



1 In the original language of |20| , the resummed prediction was based on the Catani-Trentadue quark 
and gluon "jet" functions E ?i9 defined by ]33| [m| for use in resumming logs in global two-jet event shape 
distributions. 
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supplemented by a factor containing the NGLs, 



aiC F C A -K\ 2 



R 



They identified the physical source of the additional logarithms in brackets as soft gluons 
being emitted into opposite hemispheres, with only those in one hemisphere being probed 
with the event shape pr. The additional term E px '(/7r) is generated by the terms in the 
ellipses in Eq. ( |1.4| ) that are power suppressed when m 2 2 /Q 2 ^ 1) Di rt contribute at leading 
order in Eq. fll.Tj ) because m 2 is integrated over all values. The leading double log comes 
from a light gluon jet in the m<i hemisphere recoiling against a hard qq pair in the mi 
hemisphere. Explicitly [|(], 35], 

ln 2 P fl + -- - , (1-9) 



2 In 2 2- -ln3 + 4Li 2 ( - -) + — -- 
4 I 2/ 3 6 



where the ellipses denote subleading logarithms and higher-order terms in a s . The log- 
arithm shown in Eq. Ql.9| ) is not non-global in origin and begins a series of logs that 
can be resummed by standard techniques at the leading-log level. Additional soft gluon 
emissions from the hard partons into opposite hemispheres can generate another series of 
non-global logarithms that should be included in the last term of Eq. ( |l .8[) , but they begin 
to contribute at yet higher order in a s , beyond the scope of this paper. 

Our examination of NGLs will focus on probing soft radiation in separate regions 
with different scales, by examining hemisphere masses with m 2 ,m 2 <C Q 2 using the dijet 
invariant mass distribution in Eq. . In this regime we need only consider soft emissions 
from hard q and q partons in opposite hemispheres. By simultaneously considering m 2 and 
m 2 . the two sources of large log contributions in Eq. ( |1.8[ ) can be cleanly distinguished 
at 0(a 2 ), since they arise from different parametric regimes, namely m 2 . <C tt7 2 <C Q 2 
and ml ^ m i ~ Q 2 respectively. With m\ <C m 2 <C Q 2 we have an NGL of type 2, from 



probing soft radiation in separate regions of phase space with different scales [21|. Eq. (Q 



suggests that the NGLs of type 1 and type 2 are in fact related, with the double log of pr 
in Eq. (|l.8| ) found in [20] arising from a double non-global log of h/h in the fixed order 
0(a 2 ) part of the soft function S(li,l2,fJ>) in Eq. (L4). We will use our results to verify 
this statement. 

In this paper, we calculate the full structure of the dijet soft function S{l\,l2-,p) at 
0{a 2 ) analytically, and also explore the relation to NGLs, which first arise at this order. We 
will learn that NGLs of a given observable can be understood as logs of ratios of multiple soft 
scales that are left over in the perturbative expansion of the relevant soft function even after 
a standard factorization theorem for the observable has been established. This perspective 
allows one to calculate systematically not only the leading non-global logarithms that have 
been identified in [2C, 21], but also to calculate straightforwardly the full set of other 
non-global structures that actually arise. 

We will find it simpler for many purposes to work with the Fourier transform of 
S{1\, £2, m) to position space. The cross section transforms as 

/°° 9 9 tin 

dm 2 dm 2 . (1.10) 

.oo dmfdmf; 
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In position space, the convolution in Eq. (1.4) becomes a product 



0"(yi,2/2) = cr H 2 (Q; fj.) Jnivuy) Jn(y2, m) S(Qyi,Qy 2 , ■ (1.11) 

In position space, the renormalization group evolution of each factor becomes particularly 
simple as they each renormalize multiplicatively, as opposed to convolutions in momentum 
space. In fact, RGEs in momentum space are most easily solved by first going through 



position space [36, 37 1, and there are several formalisms in use to do this to arbitrarily high 
orders as long as the anomalous dimensions are known (using derivatives on the evolution 
kernels |3gj ], converting analytically back to momentum space or even transforming 



numerically). For example, the formalism of [38| for summing logarithms in momentum 
space relies on the form of the position (or Laplace)-space soft function using derivatives 
as arguments (the "associated jet/soft functions" in p8| ). Thus the position-space soft 
function S(x±, x 2 , fi) is a crucial ingredient in resummation of logs in practically all dijet 
observables. 

We will present results for the dijet soft function to 0(a 2 s ) both in position space, 
S(x±, x 2 i A 4 ), and for the double cumulant in momentum space, <S c (^i> f- 2 i A*)- We determine 
analytically for the first time the full set of functions of x\/x 2 or If/ifi that appear in S 
and S c at this order. 2 

From our results, we learn: 

• There is a non-global double logarithm with CfCa color structure in <S c (^ij ^2> A*) 
and S(x\,x 2 , fi), which corresponds exactly to the non-global log first identified by 
Dasgupta and Salam in pO] from the single-hemisphere mass distribution. 



• There is a single non-global logarithm in S c (£f, 1% , ^) and S(x\, x 2 , /u) that appears 
when the ratio ii/£ 2 or £ 2 /£i becomes large, or when x\jx 2 or x 2 jx\ becomes large. 
Both the color structures CfCa and C^T^ny have this single log. 

• In addition, other non- logarithmic non-global structures arise in S c and S at 0(a 2 s ). 
These structures and the single log were not accounted for in previous conjectures 
about their form [fE|. 

These results not only complete our knowledge of the soft function at 0(a 2 ), together with 



known results for the hard and jet functions in Eqs. (|1.4| ) and ( 1.11 ), they make possible 
resummation to N 3 LL accuracy in these doubly-differential cross sections when £1/^2 ~ 1 
or x\jx 2 ~ 1 (when these ratios are large/small the NGLs must be resummed as well). 
In position space the complete analytic N 3 LL result can immediately be obtained from 
our results through multiplicative RG evolution, while in momentum space convolutions 
between the evolution factors and fixed order functions must still be performed. This 
becomes a more nontrivial exercise to perform in the presence of the non-logarithmic non- 
global structures found here. 



2 There is a constant term in S(xi, X2, /i), or coefficient of 5(li)8(l2) in S(h, h, ft), that we do not calcuiate, 
but which can be deduced from the recent result for the thrust soft function by ^o|. See also the note 
added to the Conclusions. 
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We arrive at the above results as follows: 



In Sec. ||, we review properties of the dijet soft function that are already known, and 
an existing conjecture for the 0(a 2 ) parts of it that are (so far) unknown. We also explore 
the relation to NGLs in detail. 

In Sec. ||, we give our final results for the calculation of the dijet soft functions in 
position space, S(xi, X2, m)> an d the double cumulant projection of the momentum space 
soft function, S c {t\ , p), to 0{ot 2 ). We organize contributions to the two-loop soft func- 
tion according to whether one or both hemispheres are populated by final-state particles. 
Thanks to renormalization group invariance and known anomalous dimensions of the two- 
loop soft function, we obtain all the new information we need from contributions with two 
particles (gluons or quark-antiquark) in the final state going into opposite hemispheres. 
These contributions encode all the non-global structure in the dijet soft function. In ad- 
dition, we will give a remarkably simple function that approximates very closely the total 
non- logarithmic non- global terms in the double cumulant S c . 

In Sec. we discuss several projections of our soft functions S(£x,£2,p) and S(x\, X2, p), 
including how we obtain analytic results for the double cumulant, S c {£\,£\, p), and the 



asymmetric thrust and heavy jet-mass event shapes of Ref. [22]. The latter include the 
standard thrust and heavy jet-mass projections as special cases, and we compare our ana- 
lytic results to EVENT2, providing a strong consistency check. 

In Sec. H we carry out another consistency check by comparing the dijet factorization 
theorem including «S c (^f , , /^) , expanded to 0(a 2 ), to the two-dimensional double cumu- 
lant distribution from EVENT2. We demonstrate that when our results are included in 
the factorization theorem in Eq. Q1.4p at 0(a 2 ) that the remaining terms are truly power 
suppressed. In this section we also study the region m^ 2 <C m'j 2 ~ Q 2 , and calculate the 
additional global double log of m^/Q 2 which arises in the double cumulant Eq. (1.5) when 



mf 2 enters the hard regime. These additional logs have the same origin as those computed 



in Ref. [35|, and are the origin of the logs of pa in the S p,c ' term of Eq. fll.SP . By projecting 



our double cumulant onto the single cumulant in pn, we reproduce the value of the double 



log given in Eq. (1.9). 



In Section ||, we conclude. 

In Appendix [AJ we record the known anomalous dimensions of the dijet soft function 
to two-loop order, which we need to assemble our result for the full soft function. 

In Appendix |b|, we provide details of the calculation of the various parts of the soft 
function given in Sec. ||[ organized by color factors and Feynman diagram topologies. We 
give results in both momentum and position space and provide details of the translation 
between the two. 
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2. The Dijet Hemisphere Soft Function 



2.1 Definitions 

The dijet hemisphere soft function is defined as 



S(£i,£ 2 



-^trV|<X s |T[Y n t^]|0) 



N, 



x s 



S[£i 



n-h \8[ £ 



(2.1 



where the trace is over colors, T denotes time ordering, = (1, 0, 0, 1) and = (1, 0, 0, —1) 
are light-cone vectors along the ±z directions, and R, L specify the ±z hemispheres which 
we refer to as right and left hemispheres. The Wilson lines Y n<n - are exponentials of soft 
gluons, 



YX =Pexp 



ds n ■ A s (ns) 



Y n = Pexp 



-w 



ds n ■ A s (ns) 



(2.2) 



where P denotes path ordering for color matrices and P denotes anti-path ordering. When 
deriving factorization theorems in SCET, soft Wilson lines are generated for each jet di- 
rection from rotating the collinear fields in the SCET Lagrangian to decouple soft gluon 
interactions with collinear quarks and gluons ||. Equation (|2.1| ) is referred to as a dijet soft 
function because Yn and Yn Wilson lines appear from the ro-jet and n-jet, and is referred 
to as a hemisphere soft function because it involves kinematic variables £±,£2 restricted to 
hemispheres. The measurement can be represented by the action of an operator, 



M(£ U £ 2 ) = 5(£ 1 -P L )5(£ 2 -P R ) , M(h,l 2 ) \X) = M{t x ,h) \X) 
where the projection operators acting on a final state give 



P L \X)= (^-k t )\X) , Pr\X) = ( 



ieR 



n-h)\X) 



and the measurement function M(£i,£ 2 ) is therefore given by 



M(4,y = s (, 



n ■ ki )S [ £ 2 



E 

ieR 



n ■ k, 



(2.3) 



(2.4) 



(2.5) 



It is possible to construct the projection operators explicitly from the energy-momentum 
tensor, which then allows the sum over states Xs in Eq. ( |2.1| ) to be removed [19]. 

A more symmetric expression for the dijet hemisphere soft function can be obtained by 
converting to the 3 representation fSIJl . Using T(Y n ) T = Y n = Pexp (ig J °° dsn ■ A s (ns)) 
and T(yJ) T = Pexp ( - ig J °° dsn ■ A s {ns)), where \ = T A A^ = -{T A ) T A^, gives 



S(£ 1 ,£ 2 ) = —{0 



ttYlY n M(li,e 2 )ylY* f 



(2.6) 



Due to the simple product structure of the position-space factorization theorem in Eq. ( 1.11| ), 
it is often easier to work with the Fourier transform of S(£i,£ 2 ), 



S(xi,x 2 ) 



d£id£ 2 e 



-il\x 



5(^1,^2 



(2.7) 
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The definitions in Eqs. ( |2.1[) and (2/7) imply that the functions S and S are symmetric 



under interchange of the two hemispheres, 

S(h,£ 2 ) = S(i 2 ,£i) , S(x u x 2 ) = S(x 2 , Xl ) . (2.8) 

Note that the function S{l\,l 2 ) has naive mass dimension —2, whereas S(x\,x 2 ) is dimen- 
sionless. 

2.2 Renormalization Group and Exponentiation Constraints 

Renormalization group invariance of the factorized hemisphere cross section implies a fac- 
torized structure for the soft function renormalization and its renormalization group evo- 
lution [E3]. Defining counterterms Zs in the MS scheme, the renormalized soft function 



which appears in the factorization theorem in Eq. (1.11) is 



S( Xl ,x 2 ,fi) = Zs l (x 1 , f i)Zs l (x 2 ^)S(x 1 ,x 2 ) haie , (2.9) 

and its inverse Fourier transform to momentum space, S(£±,£ 2 , /i), appears in the factor- 
ization theorem in Eq. (|1.4| ). Since S(xi,x 2 , fi) is dimensionless it can only be a function 
of fj,x%, fix 2 , xi/x 2 , and a s (fi), where the /i dependence is uniquely determined by the 
renormalization group. From now on we will drop the supercript "bare" for bare objects, 
which are always denoted without dependence on fj,, and we will include fi as an argument 
for renormalized functions. The renormalization group equation (RGE) for S is 

(ij-hiS(x 1 ,X2,iJ,) = js(xi,fj) + 7s(x 2 ,/u) , (2.10) 
where the anomalous dimension is given by 



js(x,n) = -Z s 1 (x,n)fi-^-Z s {x,fi) = -r cusp [a s ] ln(ze 7B x/u) + 7s[a s ] • (2.11) 



The solution of this RGE allows S at \i to be expressed in terms of S at 

S(x 1 ,X2,fi) = Us(xi,fi, fjbo)Us(x 2 , n,no)S(xi,x 2 ,fi ) , (2.12) 

where 

U s {x, fi, Mo) = e K ^'™^ (ie^xfio)^^ 3 '^ , (2.13) 

and the functions K and u are given in Appendix |A[ 

Following Ref. [^2j the \i dependence in Eq. ( |2.12[ ) can be entirely organized into the 
evolution factors by writing 

S(xi,x 2 ,n) = U s (x!, n, li Xl )U s (x 2 , n, H X2 ) eT( - Xl ' X2) , ( 2 - 14 ) 



where fx Xi = {ie lE Xi)^ 1 and the exponent T(x\,x 2 ) is independent of fj,. In Ref. \22\ 
Hoang and Kluth used non-Abelian exponentiation Q to show that the last factor 
exponentiates in this way. With fiQ set equal to fj, xi or fj, X2 respectively in the two evolution 
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factors, the second factor in Eq. ( 2.13| ) dependent on uj becomes unity, and the soft function 



Eq. (2.14) takes the simple form 

S(x u x 2 ,n) = e K ^^)+ K '^^) e n^) , (2.15) 
where Ki(xi,fi) = £T(r cusp , 7s, fi, fj, Xi ). 



Hoang and Kluth [22] also pointed out other properties satisfied by the function 
T(xi,x 2 ). In particular, it is symmetric, T(xi,x 2 ) = T(x2,x±), and there are constraints 
on the color factors in T(x\,X2) at each order in a s . For instance, in the Abelian limit 
with iif = light quarks the function T(x\,x 2 ) is one-loop exact. At two loops the only 
color structures in T{x\,X2) are CfCa and CVTrti/, where for SU(3) co i or the fundamental 
and adjoint Casimirs are Cf = 4/3 and Ca = 3, and Tr = 1/2. To order a 2 s we can write 

f( Xl ,X2) = %4 + %4 + 2^ t2( Xl /x2) , (2.16) 
where the one-loop constant t± is 

h = -C/— . (2.17) 

The arguments of the o? s can be chosen as a s (fj, Xl )a s (fj 1X2 ) or any other symmetric combi- 
nation of fi Xl 2 (e.g. ^fJ- Xl ^x 2 ) to the order we are working. The dimensionless two-loop 
function ^2(^1/^2) is unknown, and its computation is one of the main goals of this paper. 
From the symmetry x\ X2 we have t2(b) = t2(l/b). 

We will present our result for the position space soft functions by writing Eq. ( 2.15| ) 

as 

S{x 1 ,X2,n) = R(x 1 ,X2,n) + S(x 1 ,x 2 ) , (2.18) 

where all terms containing logs of {[iXi) are grouped together into R(x%, x 2 , fi) and all other 
terms are separated into S(x\, x 2 ). All quantities appearing in R or S will be expanded in 
powers of a s (fi) [in contrast to a s (fj, Xi ) in Eq. fl2.16|) ]. So R(xi,x 2 , fi) contains cross terms 
between the one-loop anomalous dimensions in and the one-loop part of T(xi,x 2 ), as 



well as terms generated by the running of a s (/j, Xi ) in Eq. ( 2.16 ). Similarly we will present 
the result for the double cumulant momentum space soft function S c {i\^2-> A 4 ) defined in 
Eq. ( |1.6| ), by writing 

SSiJ C 2,»)=K c (£i,£ C 2,ri+S c (l<i,£ C 2), (2.19) 

where 1Z C {£\, l^-, u) contains all the logarithmic plus function distributions in iij£\ or fi/i^, 
and S c {(.\^2) everything else. 3 

2.3 Numerical Results for Thrust and Heavy Jet Mass Projections 

Results for single differential event shape cross sections will depend on various projections 
of the function t2(b). In the literature results for these projections have been obtained with 



the EVENT2 generator [24, 25 1. Thrust T in e + e — > jets is only sensitive to the symmetric 



Due to some constant jU-independent terms TZ and <S C are not separately related to the inverse Fourier 



transform of R and S in Eq. (2.18), but can be obtained from these transforms by inspecting the results. 
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combination of hemisphere masses in the dijet limit, r = 1— T = (m 2 +m 2 )/Q 2 + ■ • • , where 
the ellipses denote terms that are power corrections to the dijet factorization theorem. The 
relevant projection from t2(b) for thrust is 

t 2 (l) = C F C A s[ CfCa] + C F T R n } 4 n/1 , (2.20) 

and these sj^'s contributes to the a 2 5(T) term in momentum space. Numerical determina- 
tions using EVENT2 in the literature include 

Becher & Schwartz @ : sf f ° a] = -30.0 ± 0.5 , 4" /] = 21.5 ± 0.5 , 

[CfCa] _ 0O , , -. >/] 



Hoang & Kluth [22] : s 2 FM = -29.4 ± 1.1 , a\ /J = 21.9 ± 1.5 , 



Chien k Schwartz pi : s[ CfC ' a] = -28.9 , s[ nf] = 21.7 , 



AHMSS H : s[ CfCa] = -28.33 ± 0.11 , sf /] = 21.82±g;g . (2.21) 



Recently the two thrust constants have been computed directly in Ref. [[SJ [MGL], giving 

s pc A ] = _ 28 .242 ± 0.003 , /] = 21.692 ± 0.003 . (2.22) 

For our analysis we will treat these two constants as known quantities. 4 

Another event shape of interest is the heavy jet mass, pu = max(m^, m^/Q 2 . In this 
case the only projection of t2(b) that appears at 0{ot 2 s ) is [23] 

f f t 2 (e> 9 ) = C F C A s[ C p fCa] + C F T R n f s% ] . (2.23) 

We will provide an alternate derivation of the integral moment appearing in Eq. ( 2. 23] ) 
in Sec. Numerical EVENT2 results include 

Chien & Schwartz g| : s [ 2 CfCa] = -16.6 , s£/ ] = 251 > 

AHMSS |H : sf p F ° A] = -16.79 ± 0.46 , s[ n p f] = 25.15+°;^ . (2.24) 

(No errors are quoted for the numbers in [23].) One of the benefits of our analytic results 
for t 2 (b) is that one can compute results for the constants appearing in various event 
shapes directly. In this vein in Sec. L2 we use our results to obtain analytic results for the 
combination 

f ^ [t 2 ( e *) - t 2 (l)] = (4^ - s[ CfC ^ + (4"/ 1 - s™) . (2.25) 

J 

Comparison of this combination with EVENT2 will provide a nontrivial cross-check on our 
calculations and we quote results here for easy reference 

CS@ - MGL@ : (s 2p - s 2 ) ICfCa] = 11-64, (s 2p - s 2 ) [n/] = 3.41 . 

AHMSS - MGL|| : (s 2p - s 2 ) [CfCa] = 11.55 ± 0.46 , (s 2p - s 2 ) [n/] = 3.46±g:g| . 

(2.26) 



'For further discussion on analytic results for these constants see the Note Added at the end of Sec. 
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2.4 Non-Global Logs from S(x\,X2, A*) 

In this section we make the connection between calculations of non-global logarithms and 
logarithmic terms in the function ^2(^1/^2)- The possible form of tiixxjx?) becomes sim- 
pler if we make a power expansion about the non-global limit, x\ 3> X2, which can at most 
yield logarithmic singularities. The leading singularity is related to the non-global double 



logarithm of Ref. 1 2C ] . This expansion therefore yields 



J™ (2(11/12) - sf y (|) + 4" In (|) + 4°' + (2.27) 
where the terms in ellipses are suppressed by at least one power of x\jx2- Since the limit is 

ro] 

not symmetric in x% and X2, the odd power ln(a;i/;c2) can appear. In fact the coefficient s 2 
can be extracted from the non-global double logarithm computed in Ref. [^(J. To derive 
this relation note that the double cumulant invariant mass distribution is 



2j2] 



o mV , tot , m% m. 



where only the term involving 4 is shown. In the limit Q 2 3> to^ 2 3> to^ 2 the ln 2 (m^ 2 /Q 2 ) 



term is leading, and must correspond with the double logarithm of Ref. [20| since it was 
derived in precisely this limit. The Dasgupta and Salam calculation therefore gives s 2 = 
— 2CfCatt 2 /3, or for the coefficients of the two possible color structures 



s [2 ][Cf c a] = -^ = _6.580 , 4 2][ "' ] = . (2.29) 

2.5 Hoang-Kluth Ansatz 

Hoang and Kluth [|^] argued that £2(^1/^2) should only contain powers of ln(xi/x2) based 
on the expectation that event shape distributions should contain only delta functions and 
plus distributions in the dijet limit. The log powers must then be even powers to satisfy 
the x\ O X2 symmetry and known LL results constrain the largest power to In 2 (£1/2:2) at 
0(a 2 s ). This led Ref. [§2| to make the ansatz 

tf lsatz (xi/x 2 ) = 4 21 In 2 + s 2 , (2.30) 

with two constants, S2 for thrust (given by £2(1) and hence by s 2 CfCa ^ and 4™^ from 
Eq. ( 2.20|) ), and 4^ f° r the double log coefficient. With this form the constant s 2 can be 



determined by knowing the 0(a 2 ) e)-function constants for thrust and the heavy jet mass, 
via 4 2] = 3(s 2 - s 2p )/tt 2 . The numerical results in Eq. (|]2|) give s [ 2 ][CfCa] = -3.5 ± 0.1 
and s 2 2 ^ n ^ = —1.02 ± 0.04. These results differ from those in Eq. fl2.29| ), so either the 



ansatz is incomplete or the calculation of Ref. pC| ] is missing a source of double logarithms. 
Using our full computation of t2(xi/x2) given in Sec. || we will show that it is the ansatz 
in Eq. ( |2.30| ) that is incomplete. This implies that non-logarithmic terms contribute to 
the integral in Eq. ( 2.25| ). Expanding our result for t2(xi/x2) for x\ 3> £2 we will exactly 



[21 

reproduce the result for s 2 in Eq. (|2.29[) , and derive non-zero results for the coefficients 



5 2 and s 2 
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3. Calculation of the Dijet Hemisphere Soft Function 



In this section we calculate the dijet hemisphere soft function in both position space, 
S(xi,X2, fj), and momentum space, S{l\,l%, fx), to 0{a 2 s ). To make the a s expansion we 
write 

S(xi,x 2 ,fJL) = 1 + S 1 (x 1 ,x 2 ,n) + S 2 (x 1 ,x 2 ,/j,) + 0(a 3 s ) , 
S(£ u i 2 ,n) = + StihJ^n) + S 2 (tu£ 2 ,ri + 0(al) , (3.1) 

where Si and Si are the 0{a\) terms. For nonzero ti or Xi the final results for each Si 
and Si are IR finite (see ) . We regulate the UV divergences in d = 4 — 2e dimensions 
and renormalize with MS. In momentum space the renormalization occurs through a 
convolution, while it is a product in position space: 

S(h,£ 2 ) = J dl\d£ 2 Z(£t - £[,£ 2 - £ 2 , i JL)Sgl{£' 1 ,£' 2 , n) , 
S(x 1 ,x 2 ) = Z(x 1 ,x 2 ,fi)S(xi,x 2 ,fi) . (3.2) 

The counterterms factor into single-variable pieces, 

Z(xi,x 2 , (j,) = Z(xt,n)Z(x 2 ,n) , Z(&x,£ 2 ,n) = Z{£\, fJ,)Z(£ 2 , (j,) , (3.3) 

and have similar a s expansion formulae, 

Z(x,fJ,) = 1 + Zi(xi ;/ u) + Z 2 (xi,fj) + 0(a 3 s ) , 

Z(£, fi) = 5{£) + Z x (£, //) + Z 2 (£, //) + 0(a 3 s ) . (3.4) 

The calculations can be organized by the number of cut propagators, 

C{k) = 27r5(k 2 )6(k°). (3.5) 

At 0{a s ) we have zero or one C(k) for the virtual and real emission graphs respectively. The 
0(a s ) matrix elements are straightforward to calculate and we give the results in Sec. |3.1| . 
At 0{a 2 s ) we have 0, 1, or 2 factors of C. The general structure of the calculation is 



explored in Sec. 3.2, where we show that it suffices to compute the double cut graphs since 
the remaining terms are determined entirely by the renormalization group and the constant 
s 2 . Final results for S(£i,£ 2 ,fi) and S(xi,x 2 , /i) at 0{a 2 s ) are given in Sec. |3.3| , while the 
details of the calculation are described extensively in Appendix ||. 

For easy reference we also record here the form of the hemisphere measurement func- 
tions in momentum and position space, A4^ k y(£i,£ 2 ) and M^(x±, x 2 ) for j = 0,1,2 final 
state particles with momenta k{ . It is the measurement functions that introduce the £\ , £ 2 
or x\ i2 dependence into the soft function. With no final-state particle we have a tree-level 



or purely virtual contribution and Eq. (2.5) gives 



MM(£ 1 ,£ 2 ) = 5(£ 1 )5(£ 2 ), M^(x 1 ,x 2 ) = 1. (3.6) 
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For one final-state particle of momentum k we have in momentum space 

M^ihJi) = 5(£i - k~)5{t 2 )9{k + - AT) + 5(£ 2 - k + )5{h)9{k~ - k+) 

= M[ L] (£ 1 )6(e 2 ) + 8{h)Mf\h) , (3.7) 

where k + = n ■ k and k~ = n ■ k, and in position space 

Ml\x x ,X2) = e- lXlk ~9{k + - AT) + e~ iX2k+ e{k- - k + ) 

= Mf{ Xl ) + Mf\x 2 ) . (3.8) 

Eqs. (3.7) and (3.8) show that the measurement divides into two terms, the contribution 
when the parton is in the left hemisphere plus the contribution when the parton is in the 
right hemisphere. 

For two final-state particles of momentum k\ and k 2 we have contributions from both 
particles in the left hemisphere, both in the right hemisphere, or one in each hemisphere 

Mf lM {hM - MS(4)5(4) + 5(4)MH(4) + Mg 2 (4,^) , 
M^ M (x lt X2) = A*£fi(si) +aC ] » +M [ £ R l(x 1 ,x 2 ) . (3.9) 

For two particles in the left hemisphere the functions are 

M [ k$ 2 (tx) = 5(h - k- - k~)e(k+ - k-)9(k+ - k 2 ) , 

H L 3M) = e-^+^e(k+ - K)9(k+ - k 2 ) , (3.10) 

and there is a trivial dependence on i 2 or x 2 . Similarly for two particles in the right 
hemisphere we have 

HTi^ = ^ ~ k t~ k tMK - ktW(k 2 - k+) , 

= e-^ k t+^9(k; - k+)9(k 2 - k+) , (3.11) 

with a trivial dependence on l\ or x\. For one particle in each hemisphere we have 

M[2(4,l 2 ) = 5(h - k^)5(£ 2 - k+)9(k+ - k^)9(k 2 - k+) 

+ 6(h - k 2 )5{i 2 - k+)9(k^ - k+)9(k+ -k 2 ), 

MS^L^) = e- ix ^- ix ^9(kf - k^)9(k 2 - k+) 

+ e-i^-^k+Q^- _ k+)9{k+ -k 2 ). (3.12) 

Diagrams involving J14^ LR ^ are the only ones that can simultaneously depend on t\ and £ 2 
(or x\ and x 2 ) and hence on £i/£ 2 (or x\/x 2 ). We will exploit this further in Sec. 3.2. We 
will label the different graphs by the locations of the cut partons. We call the diagrams 
where the final-state partons are all in the left hemisphere or all in the right hemisphere 
the same hemisphere terms, and we call diagrams with one parton in the left hemisphere 
and one parton in the right hemisphere the opposite hemisphere terms. 
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Finally we summarize some notation that will be used extensively below. In momentum 
space the logarithms are distributions 

„ / x rin fc ti 

£k(t)= — , 3.13 
L t J + 

while for cumulant momentum space the logarithms are 

L x = Mft/n) , L 2 = ln(^//x) . (3.14) 
In position space we denote the \x dependent logarithms as 

L\ = ln(ie yE xiij 1 ) , L 2 = ln(ie rB x 2 y) ■ (3.15) 

The non-global structure of the soft function will appear through the variables 

h a?i 1% , 

r = — , or b= — , or a = — . (3.16) 

h x 2 t\ 

3.1 0{a s ) Results 

The 0(a s ) diagrams contribute to the 0{ot%) terms, so we give them explicitly. Here 
the soft gluon can either be real or virtual. Purely virtual diagrams (no partons cut) in 
the soft function are scaleless in pure dimensional regularization, so we do not need to 
compute them. These graphs convert all 1/em's in the real emission result into l/euv' s - 
A demonstration of this which applies to our calculation here was given in Refs. [17, 47]. 
The measurement function Ai^ (^1,^2) i n Eq. ( |3.7| ) divides the calculation into two terms, 
one coming from the gluon in the left hemisphere and one coming from the gluon in the 
right hemisphere. 

The bare one-loop soft function in momentum space and position space in dimensional 
regularization are 



a s (^C F (e^fi 2 Yl r ,_,_., , , 



vr r(l-e)e 

S 1 (x l! x 2 ) = ^)C F {e-^Y T{-2e) 2e + (fe7JI )2e] (317) 

7r T(l — e) e L J 

At one-loop Z 1 (l 1 ,£ 2 , fJ,) = Z\(l\, fJ,)8(£ 2 ) + S{l\)Zi{l 2 , fi) and Z 1 (x 1 , x 2 , y) = Z 1 (xi,y) + 
Zi(x 2 ,fi) and the momentum and position space counterterms are 

^(x llM ) = ^^{-^-^}, (3.18) 



- 15 - 




Table 1: Examples of various types of loop graphs which enter at 0(a 2 s ) and their general proper- 
ties. V2 and i?2 are the 0{a 2 s ) purely virtual and purely real graphs, V\R\ are the mixed virtual-real 
graphs, and Z\ is the one-loop counterterm. For each class we indicate which of the three possible 
color structures CpjCfCaiCfTruf appear, and whether the graphs depend simultaneously on 
2^1,2 (in position space), on only a single variable at a time, or neither variable. The two-loop 
counterterm Z 2 is discussed in the text. 



which will feed into the //-dependent two-loop calculation. The remainder is the one-loop 
renormalized soft function in momentum and position space 

Sxit^M = {--A ( A S(l 2 ) - -£1 ( A 8{h) + 

vr [ n /i \nj 12 

Si(si,*2,m) = {-(2? + ID - ^} . (3.19) 

3.2 Structure of the 0(a 2 s ) Real and Virtual Terms 

The structure of the various types of 0{a 2 s ) diagrams for the calculation of S(xi,x 2 ) are 
shown in Table |]. In dimensional regularization the bare two-loop virtual graphs and 
one-loop virtual graphs with a counterterm are scaleless and need not be considered. The 
remaining terms can all be divided into single hemisphere contributions that depend on a 
function of x\ plus the same function of x 2 , and opposite hemisphere contributions that 
depend simultaneously on x\ t2 and contribute non-global terms. The bare and counterterm 
graphs can all be split into the form 

S{xi,x 2 ) = [5 same (xx) + 5 same (x 2 )] + S°w{xi,x 2 ) (3.20) 

as follows. Consider the two- loop counterterm Z 2 (xi, x 2 , //), from Eq. ( |3.3j ) 

Z 2 (xi,x 2 ,fJ-) = Z 2 (xi,n) + Z 2 (x 2 ,fi) + Z 1 (xi,fi)Z 1 (x 2 ,fj.) , (3.21) 



which is the appropriate form for Eq. ( 3.20 ), and shows that only Z\ contributes to S opp . 



Graphs that have only a single parton crossing the cut involve Al[ ,(xi,x 2 ) from Eq. 
and are part of the single hemisphere terms. This includes the mixed virtual-real graphs 
denoted by V\R\ in Table [j]. For the double cut real emission graphs the measurement 
function -M-^j k (xi,x 2 ) in Eq. ( |3.9| ) splits the result into a sum of single hemisphere and 
opposite hemisphere terms, 

R 2 ( Xl ,x 2 ) = [Rf mc (x 1 ) + RT mc (x 2 )} +R° pp (x 1 ,x 2 ) , (3.22) 
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where Rf me 



(xi) involves M [ ^ 2 (x 1: x 2 ), Rf me (x 2 ) involves M l £% 



and i?2 PP ( x i) x 2 ) 



involves A^lf^ (a?i, £2). Using results from Sec. [TT] the results can also be directly 



manipulated into the required form in Eq. ( p. 20 ). 

Next consider the color structures of various opposite hemisphere contributions. From 
Eq. ( gl2lD and Table | we see that for S opp (xi, x 2 ) both Z 2 and Z\R\ only contribute 



to Cp. Hence these terms are just part of the exponentiation of the one- loop result, and 
are not needed for our computation of the unknown terms with color structures C f Ca 



and CjrTjinf which appear in £2(^1/^2) (see Sec. 2.2). Thus the only term we need to 
consider is the double cut real emission graphs i?2 PP , ano - there are no UV divergences or 
counterterms for the color structures of interest in our calculation. 

Before proceeding to analyze results of our computation of R^ pp (xi, x 2 ), we must 
address one important issue, namely that only the sum of S opp and S sa,me terms in Eq. ( 3.2C| ) 
is infrared finite. From our computation discussed in Sec. 3.3 and Appendix [B] we find the 
infrared divergent terms in dimensional regularization are 



R 2 n>p (x 1 ,x 2 ) 



(27T) 2 
+ ••• . 



CiCA[ h + l ^ 



^ \ rp 3-27T 2 

+ C F T R n f 



%x\ix 2 [i e 



2„2 7J3 \2e 



(3.23) 



Here the ellipses denote IR finite terms and terms with the Cp color structure. The diver- 
gences in Eq. ( 3.23j ) cancel against same hemisphere contributions, and hence determine the 
IR structure in S sa,mc . In dimensional regularization the graphs contributing to <S same (xi, fi) 
give powers of (ixifiY determined by dimensional analysis and whether the graph involves 
counterterms, times functions of e. Since the same hemisphere graphs involving Z\ do not 
involve the color structures shown in Eq. ( |3.23| ), the IR divergent contributions which can- 
cel those in Eq. ( 3.23| ) can only come from i?| ame and V\R\. They are therefore uniquely 
determined to be 



\xi,x 2 ) 



" (2vr) 2 

(ixi^e 7B ) 4e + {ix 2 ^ E ) At 



l 6e 2 18e 



3 \ rp 3 

+ C F T R n f — 



2vr 2 



9e 



(3.24) 



Adding and expanding Eqs. ( |3.23| ) and ( |3.24| ) the terms involving these IR divergences 
cancel, but leave a finite remainder 



#2 PP (X1,X 2 ) + S! ame 



(xi,x 2 ) 



47T Z 3 VX2 



+ 



(3.25) 



Thus the 1/e terms in i?2 PP (^i,^2) with the form shown in Eq. Q3.23 ) can be dropped if 
we replace them by the double logarithm in Eq. ( p. 25 ). The calculation of the opposite 
hemisphere infrared finite terms in the ellipses then proceeds without further complications. 



Our result in Eq. ( 3.25| ) yields precisely the coefficient = —2ir 2 /3 of the non-global 



logarithm appearing in t 2 (xi/x 2 ) as discussed in Sec. and first computed in Ref. P0| . 
The procedure discussed here to determine the non-global logarithm was first presented at 
the SCET 2011 workshop |48fl, including the fact that it can be computed by an effective 
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theory that refactorizes scales in the soft function, and has ultraviolet divergences with the 
same coefficient as the above infrared divergences. 

Accounting for the cancellation of IR divergences, the full renormalized soft function 
can be arranged into what we will refer to as renormalized same hemisphere contributions 
that depends on \x and one of the x variables, and a renormalized opposite hemisphere 
contribution that depends on x\/x2, 

S 2 (x 1 ,x 2 ,fJ,) = [^ ame (x!,^) + ^ ame (x 2 , M )] +S° pp (x 1 /x 2 ,a,( / u)). (3.26) 

Here the double log in Eq. ( |3,25| ) is included in S?f v '{xi/x2,oi s {n)). We present our full 
result for Ek^^(x\/x2, a s (/x)) in the next section. In dimensional regularization the graphs 
contributing to S 2 ame (xx, fj,) give powers of (ixifi) 6 times functions of e, and hence only 
depend on the \x dependent logs L\ in Eq. Q3.15| ). At 0{a 2 s ) the general result is 

Sf me (xi, M ) = a 4 Lf + a 3 Ll + a 2 L? + aj, x + a . (3.27) 



Since the logarithms from infrared divergences appear in S opp all the L\ terms in Eq. ( 3.27 



are determined by the known anomalous dimensions for the soft function, through the RGE 



solution discussed in Sec. |2.2| , Hence the only thing not determined by the above general 
considerations and our opposite hemisphere calculation of 1? 2 is the CfCa and CpTj^nj 
color factor terms in the constant a®. These two terms contribute directly to the two thrust 



constants s 2 in Eq. ( 2.2CQ . We will present our final result as an analytic function of Li, 



L 2 , and b = x\/x 2 , plus s 2 CfC ' a ^ and s^™ 7 ' whose numerical values are given in Sec. \LS 
3.3 Final Dijet Hemisphere Soft Functions 



In Sec. |3.2| we showed that up to a constant the entire 0{a 2 s ) dijet hemisphere soft function 
can be computed from the opposite hemisphere double real emission graphs and the known 
0{a 2 s ) anomalous dimensions. For this calculation there are five classes of diagrams shown 
in Fig. [j]: independent emission (X) diagrams, diagrams with a single three-gluon vertex 
(T), and vacuum polarization diagrams with a gluon loop (Q), quark loop (Q), or ghost 
loop (H). Note that the C\ color structure only appears in X, and CpTjiUf only appears 
in Q, whereas the CfCa color structure appears in all classes except Q. For any graph 
whose squared matrix element is Ai, the opposite hemisphere contribution is 



/ 



^^MhMM^Cih)^) , (3.28) 



where C(k) is the cut propagator, Eq. ( |3.5D , Ai = {T, T, G,T~L,Q} and the measurement 
function J^v^j^. is given in Eq. ( 3.12j ) for both position space and momentum space. We 



present details of their calculation using Feynman gauge in both position and momentum 
space in Appendix |b|. The results appear in the following equations: 



X: EqB.( fT2g , g32j , p27| fc ggg ) , T : Eqs.( p^2] gj3p , Q : Eqs.( |H40| , p^ ) 



U : Eqs.( [M40|jFl43D , Q : Eqs. (p^7]p3oD . (3.29) 
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a) b) c) 




Figure 1: 0{a 2 s ) opposite hemisphere diagrams. The endpoints of the gluons can be attached to 
the points on the Wilson lines labeled by a 'x' in any order. Figure (a) gives the I diagrams, (b) 
and (c) give the T diagrams, (d) gives the Q diagrams, (e) the H diagrams with ghosts, and (f) the 
Q diagrams with massless quarks. 



We have explicitly checked that our final result is unchanged if the gluon propagators 
in Fig. [I] are taken in a general covariant gauge. The gauge parameter cancellation occurs 
individually for the T, Q + H, and Q terms (and provides a non-trivial cross check on the 
relative overall signs of Q and %). 

Next we present final results for the renormalized soft function that includes contribu- 
tions from both the same hemisphere and opposite hemisphere terms, using the approach 



described in Sec. 3.2. We first discuss position space and then the double cumulant distri- 



bution in momentum space. Eqs. ( ^30|) and (|3.36|) are the main results of this paper. 



3.3.1 Result in Position Space 

In position space we find 



S(xi,X2,n) 
where 



1 



a s (n)C F 2 , ~ t _ _ A , a 2 s {n) 



Air 



-it + R(xx,X2,/j) + 



4vr 2 



7T 



c F — + -h[ — 



t2 ( x -i) = -c F c A 2 ^- i^m 



X 2 < 

2 |q ( X 1 /X2+X2/X 1] ( CfCa ^ 2 -^C3 + CpTRnf ^l 



(3.30) 



(3.31) 



+ 2C F T R n f 
+ 2C f Ca 



9 

2Fq(1) 



F N (^) + F N {^-) - 2F N {1)] + C F C A sf FCA] + C F T R 



m f s 



f s 2 



determining the non-global function appearing in Eq. ( 2.1q ). Here and throughout this 
paper x\ and X2 have a small imaginary components, and should be regarded as xi — i0 + 



-19- 



" 15 

+ 5 

X 

£ 

CD c 





[ 1 1 J 1 1 I 1 1 I 1 1 I 1 1 [ 1 1 I 1 

\ — Ca i 

\ nf 


. \ 

■ 1 1 1 1 1 1 1 1 1 1 





-40 -20 20 40 60 80 100 
Xi/x 2 

Figure 2: The non-logarithmic non-global function in position space. We show the functions 



Re[F(fe) + F(l/b)), where b = x\/x2, that appear in the position-space soft function Eq. (3.30) for 
the color structures CfCa (solid red) and CfThu/ (solid blue). 



and x 2 — i0 + . Values for the constants are given in Sec. 2.3. The functions -Fq.tv ar e 



Fqib) 



F N (b) 



b In 2 b 3-2tt 2 



2 In 6 

3(6-1) ~3(6-l) 2 9~ 

vr 4 ln& b In 2 b 3-llvr 2 + 18C 3 / 1 
~36~3(&-l) + 6{b-l) 2 + 18 n V + b 

vr 2 T r _ ,, n2 22 



ln(b+-)+- ln 2 b ln(l-6) + -ln& Li 2 (6)-4Li 3 (b), 



11 9 , \ In b 
In b ln(l-6 ) + 

6 v ; 24 



Li 2 (l-6) + [Li 2 (l-6)] z -^ln6 Li 2 (6) + 21n6Li 3 (l-6) + 11 Li 3 (6), (3.32) 

and the formula for the //-dependent terms, R(x\, x 2 , fi), will be given below. The result in 
Eq. ( [3.30 ) is written so that the double and single logarithmic singularities for b = X\/x% — > 
and b — > oo are separated out in the first two lines. The more complicated structures 
in -Fq(6) + Fq(1/6) and Fjy(b) + F/v(l/6) are bounded on the real b axis and have real 
parts going to zero for b — > and b — > oo. For negative Re b, the functions also have an 
imaginary part. The real parts are illustrated in Fig. 0. 



Taking the x\ x 2 limit discussed in Eq. ( 2.27] ) yields terms ln l (xi/x 2 ) at 0\ 



and our result in Eq. ( |3.3C ) determines the coefficients 

.2 



s [2] ~ 
s [0] ~ 
where 



~c F c A , 4 11 



9 



(3.33) 
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-af In 2-C F C A 2F N {\) -C F T R n f 2F Q {\) + C F — +C F C A 4 



2 /4tt 2 2\ 
2F Q( l) = - + (---)ln2-8 C3 , 

2F N {1) 



[CfCa] 



C F T R n f s 



1 

3 



vr 4 I IItt 2 

-+(-- — + 2&)]n2 + 22&. 



f*2 > 

(3.34a) 
(3.34b) 
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Figure 3: Non-global logarithms and non-global non-logarithms in position space as a function of 
Xi/%2 (left panel) and in cumulant momentum space as a function (right panel). We plot 

separately the log 2 (dotted), log (dashed), and non-log (solid) parts of the various functions that 
appear in the soft functions S{x\,X2, jtt) and S c {i\, t\, fi) (ignoring the fi dependent terms in R and 
7Z C ). The single log curves have been plotted with the opposite sign. 



In the left panel of Fig. ||] we plot the double log (dotted), single log (dashed), and non-log 
(solid) 0{a 2 ) non-global contributions to S{x\, X2, //) as a function of x\/x2- Results are 
shown separately for the CfCa and CpTpjij color structures with overall signs for the plot 
chosen so that the former are negative and latter positive, which involves plotting the single 
NGL terms oc ln(6 + 1/6) with the opposite sign. The solid lines are the non-logarithmic 



functions [F N>Q (b) + F N>Q (l/b) - 2F N>Q (1) + 



1 [C F C A ,n f h 



2 S 2 



the dashed lines are the single 



NGL terms, and the dotted black line is the double NGL oc In b (nonzero only for CfCa), 
where b = x\jx2- For values in the vicinity of x\ = %2 the contributions from all sources 
are comparable in size and the non-analytic non-log dependences are important. As we 
already remarked, the F-dependent terms exhibit a smooth bounded behavior. 

Expanding Eq. ( |2.15 ) and inserting the known anomalous dimensions the L\ t 2 = 
\si{ixi t 2& lB t-t) dependent part of the soft function, R(x%,X2, (J,), is 



R(x 1 ,x 2 ,n) 



*sWC F (J 2 + zi) + ^ic F (8Lf + 8Ll + 16ZfZl) 



+ 



+ 



7T 



16vr 2 



32 



+ ( -—C F C A + ^C F T R n f ) (L{ + L 



4tt 2 C 2 f - 4C F C A 



9 



67 7T Z 



, ^0 (~t rp 



. 808 22vr 2 \ ^ m /224 8 2 

( r( A i - -jf ~ — + 28C3 j + C F T R n f I — + - tt 2 



(3.35) 



3.3.2 Result in Momentum Space 

Our complete result for the bare two-loop opposite hemisphere piece of the momentum 
space soft function S{i\,i2) is given in Eq. (|B.55[) . In Sec. 4T we describe how this result 
can be used to obtain the double cumulant S c (£^, fx). The final result for the renormalized 
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momentum space double cumulant distribution is 



4?r 



Ot s (n) 2 ( 7T 



+ 



11^-3-18(3 6-4^ 
CfCa ^ h C F T R nf 



9 



9 



In 



C F T R n f 



An 2 

2 



'2p 



(3.36) 



where the logarithmic dependence is isolated on the first two lines and the remainder 
depends on the constants s^^'™^, /i-dependent terms 1Z c (l\, l\, fi) to be discussed below, 
and the functions 

/Q(«) =(^Y~ 3(^1)) ln (°) " \ ln (°) Li2 (" a ) + 4Li 3("«) " " 2^ 2 ) In (a + ~ 



/at (a) = -4Li 4 

+ Li 2 
1 

+ 24 



1 



a + 1 
1 

a + 1 
22 In 



HLi 3 (-al + 2Lij 



1 



a + 1 



In 



vr 2 - In 2 (a + 1) - -ln(a)ln 



1 



(a + 1) 2 
a 



L(a + 1) 



11 1 / N 

+ y ln(a) 



(a+1) 2 



6 In 1 + - ln(l + a) + it 1 } ln 2 (a 



(a - l)ln(a) 
6(a + l5 



5-7T 2 / 1\ , . 

+ l2- ln ( 1+ a) ln(1 + a) 



117T 4 (11^-3-18(3) , / , 1 

In a H — 

180 18 V a 



(3.37) 



The combinations /g(a)+/Q(l/a) and //v(a)+/jv(l/ a ) appearing in Eq. ( |3.36D are bounded 
functions of a, vanishing as a — )• or a — )• 00. 

Taking the l\ 3> ^2 limit of Eq. ( |3.36j ) we see that the non-global logs in momentum 
space have the same coefficients as those in position space. We find 



8^2 



4 2 hn 2 m+4 i jinm+4 u j + 



'2c 



'2c 



'2c 



(3.38) 



where s 2 2 J = and 4c = s 2 ±J from Eq. ( 3.33 ). The constant s 2 c> meanwhile, is given by 



7T 



4c = -4c ln2-C F C A 2/^(l)-C F T i? n / 2/Q(l) + C F — + OC^, + C F T fi n /S2p 



where 



/] 

(3.39) 



2/ Q (l) = -6C 3 --(3-2vr 2 )ln2. 



(3.40) 



2/jv(1) = -8Li 4 - +C 3 hr-51n2 + 



33 



ln2-ln 4 2 2tt 4 vr 



+ — hi z 2 



11 



ln2 . 



3 45 3 V 3 

In the right panel of Fig. ^ we plot the double log (dotted), single log (dashed), and non-log 
(solid) non-global contributions to S c (£^, £ 2 , jj). The overall signs for the plot are chosen 
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Figure 4: The non-logarithmic non-global function in momentum space. We show the functions 



/(a) + /(1/a), where a = t%/t\, that appear in the momentum-space double cumulant Eq. (3.36) 



for the color structures CfCa (solid red) and CpTjiUf (solid blue). The simple approximate forms 



Eq. (3.41) of these functions are the dotted lines. 



so that the CpCa color structures are negative and CpT^rif are positive, which involves 



plotting the single NGL oc ln(a + 1/a), where a 



?-f , with the opposite sign. The solid 



lines are the non-logarithmic functions [/jv,Q(a) + /jv,q(1/o) — 2/jv,q(1) + \s 



i [C F C A ,n f h 
2 b 2p J- 

Comparing to the left and right panels of Fig. [|| the log terms are identical, but the non-log 
terms differ. 

The symmetric sum of complicated functions in Eq. ( |3.37| ) that appears in Eq. ( 3.36|) 
is well approximated for all a by a simple function, taking the form 



4a 



>M a ) + fQ^f (1/a) n 2/ Q ,jv(l) - + . 



(3.41) 



The comparison between the exact forms in Eq. ( p. 37 ) and the approximate forms in 
Eq. (3.41) is shown in Fig. |], displaying remarkable agreement. 

and s^J^ in Eq. (|3.36|) are given by the thrust constants 



The constant terms s 



JCfCa] 



[C F C A ] 
2p 



and , plus analytic constants. In Sec. 4.3.2] we calculate the differences s% — sjf 



analytically, and find 



»2p 



SC F C A ] 
5 2p 



2 3 45 3 3 V 2 



(3.42) 
HC 3 -14C 3 ln2. 



Thus only the constants s 



[C F C A ,n f ] 



remain as unknowns in both of our final results Eqs. ( |3.30 ) 



and ( |3.36| ), and for them we use the numerical values quoted in Eq. (2.22) from Ref. |fi"o|j . 
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The ln(£i /fj,) and ln(^|///) dependent part 1Z c (£i, £$, A*)> obtainable by taking the in- 
verse Fourier transform of Eq. ( 3.35| ), is given by 



n c (£i,£ c 2 ,fi) 



32 



a s (/x)C. 



7T 



jsC 2 (L? + 4) + leC^L 2 . + (^C F C A 



-C F T R n } ) {L\ + Ll) + 



207T- 



(Li + ^) 



+ 



64C3C | + C F C A ( f - - 28C,) - c F r sn/ (f - ^ 



(i>i + i> 2 ; 



(3.43) 



where = ln n (^f///), and L% = ln"^/ A 4 )- Eq. ( |3.36| ) will be used in the double cumulant 
comparison to EVENT2 data done in Sec. [|. 



4. Projection onto Other Observables 

Having obtained the explicit analytic result for the dijet soft function, we can make use of 
it in several ways. In Sec. |4.1| we give further details about the projection of the momentum 
space result used to obtain the renormalized double cumulant, 6 ' c {£\, £%■> m) m Eq. ( 3.36p , 



since the procedure is useful for other projections. Additionally, in Sec. 4.2 we project 
the distribution onto two classes of e + e~ event shapes r Q and p a originally defined in J22[| , 
which we call asymmetric thrust 5 and asymmetric heavy jet mass, respectively. Asymmetric 
thrust can be defined from the hemisphere masses and m| as 



amf + 



(4.1) 



a 1 + a Q 2 

where a > is a dimensionless parameter. In the dijet limit r a corresponds to the original 
thrust for a = 1. The asymmetric heavy jet mass is 

2 



(4.2) 



and corresponds to the heavy jet mass for a = 1. The da/dT a and da/dp a distributions 
contain non- global logarithms in a, and we will show that the full non-global structure of 
their cumulants bears close relation to the non-global structure in the position-space dijet 
cross section and the double cumulant momentum-space cross section. 

As discussed in Sec. |2.3| , the projection of the non-global structure onto the two-loop 
constants needed for the heavy jet mass and thrust is of interest. These observables have 
recently been used in fits of a s from e + e~ collider data, placing importance on the knowl- 



edge of these constants [16, 18, 23, 5Cj. In Sec. 4.3 we project our non-global results onto 
heavy jet mass and thrust and determine the difference between the cumulants of the heavy 
jet mass and thrust distributions, Eq. ( |2.25| ). We compare our analytic result, Eq. (|4.29| ), 
to a numerical projection from our position-space result and to a numerical extraction from 



EVENT2, Eq. ( p. 26 ), both of which provide a nontrivial cross-check of our calculations. 

5 This is to be distinguished from the hadronic event shape ta for pp collisions called asymmetric thrust 
in M. 
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4.1 The Double Cumulant of the Non-global Terms 

To determine the double cumulant of the non-global terms, we compute the opposite hemi- 
sphere cumulant 

s c opp (^V 2 c ) = / 1 Mi / 2 de 2 s opp (h,e 2 ), (4.3) 



The function S opp (£i,£ 2 ) only has support for £\,£ 2 > 0. The relation of this double 
cumulant soft function to the position-space soft function is 

sraiJi) = - r PP<-*ul* + , (4.4) 

J_oo 27T 2lT {XlX 2 ) 

where x% = %% — i0 + in the above equation. Comparing this result to Eq. (4.14) below 
with a = 1 we see that the diagonal projection S ovv {£\,£\) is equal to the heavy jet mass 
cumulant with £\ = Qph, which explains the appearance of the heavy jet mass constants 
JCF c A,nj] q ^ final result in Eq. (|3.36| ). 

For the analytic computation of the dependence on £\j£ 2 of the opposite hemisphere 
terms <S° pp (^f, £ 2 ) it is more convenient to start with the bare momentum-space result. We 
use the general form of the opposite hemisphere results from Eqs. ( |B.13j ) and ( |B.17| ): 

S°ppM(^ 2 ) = AC i ^ S - 1 ~ 2e [F + F 1 (r)] , (4.5) 

where the variables 

s = , r = l f , (4.6) 

the color factor Cj = CfCa or CpT^nf, F and i*i(r) depend on the class of diagram (see 
Fig. 0) contributing to 5opp and with \i in MS the prefactor 

To evaluate the double cumulant we change variables to s and r: 

S c opp (^ 2 c ) = / - / dsS°™(£ 1 ,£ 2 ) + / - / dsS°^(h,£ 2 ). (4.8) 

Jo 2r Jo Jq/t? 2r Jo 

Plugging in the form in Eq. ( |4.5| ), we can perform the s integral. Then expanding to collect 
terms, we find 



WM) - C [±A Fo (£f - (£)' 1 [ ± Fllr) ,4,, 



2 V2vr/ J n r Wf / 2 V2tt/ 



The final two terms contain only finite //-independent terms which contribute to S. The 
first two terms contain //-dependent terms, some of which cancel with the same hemisphere 
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terms in the same way as explained in Sec. 3J2. The only term that does not cancel, but 
combines in a non-trivial way, is the double logarithmic term. Putting these terms together, 
the cumulant of the non-global terms S is 



(4.10) 



+ 



1 WW dr, ( l 



In (r^)if ) (r) + 



1 dr ( i 

— In 



Fi 0) (r)+Sa 



where Sq is the constant term in the non-global terms and is given by 



~ _(o) F< 



(2) 



S " = "24 F » ' + 



1 f°° dr f 1)t . 



(4.11) 



and is the coefficient of e n in the expansion of -Fq,i- We can perform these integrals 



analytically for the total momentum space result from all graphs, given by Eq. ( B.55| ). 



Note that our calculation of 5° pp ^ CfCa ^ gives 



F, 



(o) 



2tt2 
3 



(4.12) 



The final result for the double cumulant with both opposite and same hemisphere terms is 
given above in Eq. ( [3,36 ), 



4.2 Asymmetric Heavy Jet Mass and Asymmetric Thrust 

In this section we give the projection from both the position and momentum space opposite 
hemisphere terms onto asymmetric heavy jet mass and asymmetric thrust. We will use 
the a = 1 results of these projections to extract the difference in the constant terms 
of heavy jet mass and thrust. Projecting both forms of the hemisphere soft function 
result, position and momentum space, onto the asymmetric heavy jet mass and asymmetric 
thrust cumulants provides a nontrivial check of the consistency of our results and uncovers 
interesting relations for the hemisphere soft function. We present the relevant projections 
starting from the hemisphere soft function in both spaces. 

The asymmetric thrust cumulant is defined from the momentum space and position 



space hemisphere mass soft function by the projections [22] 



dr n 



dmi drrin S 



2 ami + m?. 



1 + a 



Q 2 



d 2 a 



dmfdm2 



2 ' 



dy 

— exp 
2vr P 



iQ 2 r c a y- 



y-i0+ 



(4.13) 



Similarly, the asymmetric heavy jet mass distribution is defined from the momentum space 
and position space hemisphere mass distributions by p2| 



dp, 

dyi dy 2 
2vr 2vr 



dm\dm\ 5 



max^amj, m 



exp 



iQ 2 p c a (yi + y2 



1 + Q 

(1+a) 



Q 2 



d 2 a 



dm^dm^ ' 



&(ayi,y2) 



(yi - i0+)(y 2 - i0+) 



(4.14) 



- 26 - 



+iOO 



S-Jt -2 -1 



z>o 







-it 



-ICO 
— 



z<o 







■Jt -2 -1 



-lOO 
— » — 



3 
2 
1 

-1 
-2 
-3 



+lOO 







i 2 +jt: 



Figure 5: Original integration contours for 9 in the complex plane after the change of variable in 
Eq. Q. 



4.2.1 Projection of Non-global Terms in Position Space 

Let's consider the contribution to the two asymmetric event shapes from the non-global 
term. In position space we can determine this contribution by inserting ^{yi,y 2 ) = 
a 2 s /(8TT 2 )t 2 ( yi /y 2 ) into Eqs. (fDJ|) and flOg) . 
For asymmetric thrust this gives 



( 



o 

a 



* exp (i Q^ s) -*»« 
2vr y 



i0+ 



(4.15) 



This result tells us that the non-global dependence on a of asymmetric thrust, given by 
t2(a), is equivalent to the non-global dependence on x\jx 2 in the position space hemisphere 
soft function, given by t 2 (xi/x2)- We will see a similar correspondence for asymmetric 
heavy jet mass in the momentum space projection. 

For asymmetric heavy jet mass the fact that t 2 only depends on the ratio of positions 
motivates making the following change of variables in Eq. Q4.14| ): 

i0 + .„j_ dy\dy 2 idO dz 



id = In 



yi 



yi + yi - tif 



(4.16) 



y 2 -i0+' " a± 1 "* " (lft-»0 + )(l&-»0 + ) (z-iO+)' 

With the new variables z is integrated along the real axis, — oo < z < oo, and 9 is integrated 
along a contour in the complex plane shown in Fig. || which depends on the sign of z. For 
either sign of z this contour can be deformed to the real axis without encountering cuts or 
poles, so that the limits become —it < 9 < tt. This gives 

„2 /■+* d Q ,+OC ^ 



^PaKPa) g7r 2 / 27r / 27r 



exp (iQ 2 p c a z) t 2 (ae ie ) 



(4.17) 



The projection again just gives a ^-function of p^,, but this time a more complicated 
projection from t 2 is needed to get the coefficient. 
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Note that the calculation used to obtain Eqs. ( |4.15 ) and ( 4,17 ) is also valid for terms 
in the cross section at higher orders in a s as long as they are only functions of x\/x2- 



4.2.2 Projection of Non-global Terms in Momentum Space 

For both asymmetric heavy jet mass and asymmetric thrust, we consider the projection 
for the bare opposite hemisphere terms and remove the //-dependent pieces to obtain the 
projection for the non-global terms. To perform the projection, we change variables from 
i\ii-2 to s and r in Eq. (|4.6|). For asymmetric heavy jet mass, the cumulant for the opposite 
hemisphere terms is 

Sc°Z P (P«) = £" d Pa J dhd£ 2 5 (p a - ^— max |)) S°™(h,£ 2 , p) . (4.18) 

Using the general form in Eq. ( [4.5[ ), the cumulant in terms of integrals over s and r is 

2 2 p- 2 ((l+ Q )/2a) 2 r 



( f a tir f 

ScZ pb] (P c a ) = d(P c a ){J o Yr I 



u,/ ; , , ^ 4e„— 1— 2e 



+ J — J dsAC i ^ e s- 1 - 2e (F + F 1 (r))j. (4.19) 

From this formula we see the p a cumulant is just the momentum-space double cumulant 
Eq. (fOl) with £± = Qp c a (\ + a) /(2a) and £| = Qpafi + ol)/2. Performing the s integral in 
Eq. ( 4.18| ), expanding up to O(e ), and canceling 1/e and //-dependent terms against same 
hemisphere pieces as in Sec. and Eq. ( 4.10 ), we find 

Sc™ = 0( P l) (g) 2 Q| - \F®> In 2 a + S (4.20) 

1 /• 1/Q dr ln(ra)F ( 0)(r)+ 1 r± Hr/a)F (v) {r) \ 

2 Jo r 2 J r J 

From this result we find the non-global structure in the last two terms is equivalent to the 
corresponding non-global structure in the double cumulant, Eq. ( 4.10| ), with the replace- 
ment l%/l\ — > a. The constant So is the same as in Eq. ( 4.1 1| ) . 

We can go through the same procedure for asymmetric thrust, whose cumulant is 

Scrl p (r c a ) = jf" dr a j dhd£ 2 5(r a - ^-^1±^\ S^^^p) . (4.21) 
In terms of integrals over s and r, 

ScT«) = l Yrl dsS ° pp ^^ri> ( 4 - 22 ) 



where 

(Qr^(l + af 



Aa[ar + l/(ar) + 2] 



(4.23) 
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Evaluating the s integral, expanding to finite terms, and canceling [i- and e-dependent 
terms against same-hemisphere contributions as in Sec. [T|and Eqs. flOcD and ( p0|) , 



S™ « ( T 3 =*( T o) (^) ! C,{ 1 In 2 a - i to (2 + or + ^) ff (r) + So, } . 

(4.24) 

where So T = 5o — (vr 2 /6)F (0) . Although the non-global structure in this result is less 
obvious, we know from the position space result that the non- global structure in a for r a 
is functionally equivalent to the non-global structure in t 2 (xi/x 2 ), the position space form 
of the non-global terms in the soft function. 

To conclude, we have shown that the non-global structure in the hemisphere soft 
function in position space and momentum space is reproduced by the asymmetric thrust 
and asymmetric heavy jet mass observables respectively. 

4.3 Heavy Jet Mass and Thrust 

For heavy jet mass p# and thrust T = 1 — r we will compute the difference of the cumulants 
of the non-global terms at 0{a 2 s ). This corresponds to the projection of t 2 {x\,x 2 ) in 
position space, Eq. ( |2.16 ). From Eqs. ( 4.1 5| ) and ( 4.17| ) with a = 1 this simply gives 



constant terms in the cumulants, which are S(pn) and 5(t) terms in the corresponding 
distributions. Since we will take the difference of heavy jet mass and thrust we only need 
opposite hemisphere terms. We examine the projections in both momentum space and 
position space. The results are the same, providing an internal consistency check on our 
position and momentum space results, and more importantly the analytic results for the 
projection match the numerical extraction of these constants from EVENT2, which were 
discussed in Sec. ^j^. We call the heavy jet mass cumulant and S r the thrust cumulant, 
and add a superscript NG to refer to keeping only the terms in the cumulants induced by 
the non- global structure at 0(a 2 ) 



4.3.1 Projection From Position Space 



Taking the a = 1 result for the asymmetric heavy jet mass projection from Eq. ( 4.17| ) 



to find the projection for the non-global terms and plugging in the result Eq. ( p. 31 ) for 

t 2 {x\/x 2 ), we find 

8tt Jo vr (425) 
= ^ 0(P) [c F C A sf; CA] + C F T R n f s [ % ] ] , 

where we usedt 2 (e ie ) = t 2 (e- ie ). Similarly, we can take the a = 1 result for the asymmetric 
thrust projection from Eq. ([4.15|) : 



a 2 „, _ a 2 



^r G (.r) = ^e(T)t 2 (l) = ^d(r) C F C A s 2 CFCAl +C F T R n f s 2 nfi . (4.26) 



[C F C A ] 
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Therefore the difference in cumulants is 



5^ 



C F C A {s 2p - s 2 ) [CfCa] + C F T R n f (s 2p - s 2 ) 



(4.27) 



V] 



which provides an alternate derivation of Eq. ( |2.25| ). We will use the momentum space 
projection to determine the analytic value for the difference in cumulants, as the integrals 
are more straightforward. 

4.3.2 Projection from Momentum Space 

In momentum space, we can take the a = 1 results for asymmetric heavy jet mass and 



asymmetric thrust, given in Eqs. ( 4.20 ) and ( 4.24 ) respectively. When we take the differ- 
ence, what remains is the difference in constants for the two distributions: 

.2 ,1 dr 



iNG 
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^NG 



(P) 



2?r/ 



) 0{p)Ci 



■Fo + 2 



ln(l + r)Fi(r) 



(4.28) 



a 



^e(p) C F C A (s 2p - s 2 )^ c ^ + C F T R n f {s 2p - s 2 ) [n/J 



where the s 2p and s 2 constants are defined in Eqs. fl2.20|) and ( p.23|) . Note that — \Fq gives 
the coefficient of the non-global double log term in the hemisphere soft function. Using the 
results in Appendix M we find 
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llC 3 -14C 3 ln2 



(4.29) 



This yields the analytic result quoted above in Eq. ( 3.42] ), and agrees with a numerical eval- 



uation of the projection in Eq. (4.27) using our analytic position space result. Numerically 



Eq. (|4.29| ) gives 



{s 2p - s 2 ) [CfCa] = 11.6352, (s 2p - s 2 ) [n/] = 3.4749 



(4.30) 



which can be compared to the EVENT2 extractions in Eq. ( p. 26 ). There is excellent 
agreement between our analytic results and the EVENT2 extractions, providing a strong 
consistency check on our calculation. 

The constant s 2p — s 2 has contributions from both the non-global logarithms and the 
remaining non-global terms. For each color structure, the contributions from each set of 
terms is 

C F C A C F T R nf 

double log: 43.293 - /in ^ 

, , (4.31) 
single log: 39.586 -15.791 v 1 

non-log: -59.609 22.741 . 

We see that the non-global logs do not dominate the contribution, and the full set of 
non-global terms is needed to accurately determine the constant. 
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5. Two Dimensional Comparison with EVENT2 



The Monte Carlo EVENT2 p4| , |25| allows us to study the structure of the hemisphere jet 
mass distribution. EVENT2 contains the matrix elements necessary to compute dijet event 
shapes at O(o?), and it allows us to numerically compute the non-global contributions 
in the hemisphere mass distribution. We use the program to numerically compute the 
distribution of m\ making the cut m\ < m^ 2 : 

da 'mf)= r f dml-^. (5.1) 



dm 2 , J dm 2 dm 



In the region m\ <C m^ <C Q , the dijet factorization theorem we use here applies, and 
we can compare the results of EVENT2 in this regime to our calculation. In Sec. 5.1 we 
compare the soft non-global logarithms from EVENT2 with our calculation and show there 
is excellent agreement. 

When ^ m< \ ~ Q 2 ' our factorization theorem does not apply and new contri- 
butions arise. In Sec. [5^ we consider the logarithmic structure in this hard region and 
compare the leading double log to EVENT2. 

5.1 Soft Non-Global Logarithms 

In this section we perform a direct comparison of our hemisphere soft function with the 
Monte Carlo program EVENT2. The distribution in Eq. ( |5.1| ) is related to the double 
cumulant in m\ and m 2 , S(mf 2 , mSf )■ For a bin in m\ between m\ min and m| max , the 
value of the distribution in the bin will be 

S(mf,m| max ) -S(mf,m| min ) 



„ / 2, max „ d^(J / 1 , 9 / 2 ' min ,0 (1 (7 

2 A >m 2 



dmi / d ™ 2 An 2 At; 2 / / Ar7 2 r7m 2 ' (5 ' 2) 

o jo urri]^ am 2 Jo Jo am 1 ara 2 



This double cumulant is for the full distribution, meaning it is a sum of the global terms in 
the distribution, the non-global terms that we have determined, and terms that are power 
suppressed by m 2 /Q 2 « 1. To test our analytic results we will make use of the double 
cumulant for the non-global terms given in Eq. ( 3.36| ). The region we are interested in to 



avoid power corrections and obtain sensitivity to non-global logs is 

m\ < mf < Q 2 . (5.3) 

As m\ — > Q, we become inclusive over a region of phase space, m\ <C m^ 2 ~ Q, where 
additional logs of m\jQ that were power suppressed for m\ <C Q begin to contribute. These 
additional logs are non-singular in mi, e.g., m,i In(m2/Q). As we vary m^ 2 in the program, 
we observe the effects of such terms, and we can turn them off by choosing mf 2 <C Q 2 ■ 

In our comparison to EVENT2 we focus on the CfCa and CpT^rif color structures, 
as the form of the Cp color structure is constrained by exponentiation. To compare our 
calculation with the soft non-global logarithms in the double hemisphere mass distribution, 
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we subtract the full dependence in these color structures for the double hemisphere mass 
distribution from EVENT2. 

The level of numeric accuracy in the doubly differential distribution required to resolve 
the contribution of the non-global terms is significant. To achieve this accuracy for the 
C F C A terms we ran EVENT2 with 1.82 x 10 12 events. Additionally for the C F C A color 
structure, we chose the parameters CUTOFF = 10 -15 and NPOW1 = NPOW2 = 3. 
The choice of CUTOFF allows us to reliably sample m\ down to log 10 tti 2 ./*? 2 = 10~ 7 
and the choices of NPOW1 and NPOW2 significantly reduce the statistical uncertainties 
compared to the default choices NPOW1 = NPOW2 = 2. For the CpT^nf color structure, 
we found that taking NPOW1 = NPOW2 = 4 was necessary to obtain reliable results, 
and with this parameter choice we ran 0.253 x 10 12 events. This choice of parameters, 
though, introduces large systematic and statistical errors into the large m\ region of the 
distribution. It appears to be more delicate to compare the results of EVENT2 for this 
color structure with our calculation 6 . 

The full set of terms contributing to the distribution can be divided up into groups. 
In the limit we are working, the cumulant terms can be divided up as 

S = £ global + S ln2 + S ln + S NG + S px - . (5.4) 

The global singular terms determined by RG evolution in the dijet hemisphere factorization 
theorem, are in s global . The non-global double and single log terms are in S ln2 and S ln , 
and the remaining non-global non-log terms are contained in E NG . Here S ln and S NG are 
determined by our calculation. Additionally, there are power corrections to the dijet limit 
that sit in S px ' and contribute in EVENT2, but which are not known analytically. 

For both of the color structures CpC A and CpT^rif we will show the distribution in 
m\ from EVENT2 with a cut on m 2 , da E V2/ 'dm^m^ 2 ) , with various terms subtracted. 
We choose the cut in m% to be 

mf/Q 2 = 9.2 x 10" 4 , log 10 (mf /Q 2 ) = -3.0 . (5.5) 

Specifically, in Figs. |6| and we plot the distributions 

(5.6) 



(top-left) 


dm 2 (mi } ' 






(top-right) 


dm 2 (mi } 


do-global (im c2\ 

dm 2 {mih 




(bottom-left) 


dm 2 (mi } 


""global ( c2\ 

dm 2 {mi) 


da ln 2 
dm 2 . 




dm 2 (mi } 




dain ( 
dm?. 


(bottom-right) 


da E y 2 , c2 s 
dm 2 (mi } 


^global ,_c2\ 

dm 2 (mi) 


dm 2 . 



mf) for C F C A 



(mf ) for C F T R n f 



mf)-p^(mf)- d -p^(mf) 
dm% dm% 



6 The four-particle phase space generator in EVENT2 is optimized to sample the singular regions of phase 
space for a qqgg final state J5l]]. This may be a reason why the CfTjjti/ results, which have contributions 
from matrix elements with a qqqq final state, are less consistent. 
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Figure 6: The Q{a 2 s ), CfCa distribution from EVENT2 with various sets of terms removed. The 
upper left plot shows the distribution directly from EVENT2. The upper right plot has the known 
global terms subtracted, and the lower left plot has both the global terms and the known non-global 
double log in the cumulant subtracted. The lower right plot has the global terms and the non-global 
terms from our calculation subtracted. The red line marks the value of m° used. Note the changing 
vertical scale between the plots. 



In Fig. ^ we show the distributions in Eq. (|5.6| ) for the CfCa color structure. Note that a 
non-zero flat region in the distribution indicates the presence of a single log of m^ 2 , while 
a line with non-zero slope indicates the presence of a double log of m^f. We can clearly 
see that after we subtract the global terms a double log remains. Once we subtract this 
double log coming from S ln , the remaining distribution still has a single log of m^ 2 in S ln . 
Once we remove it and the remaining non- logarithmic non- global terms that are in £ NG , 
the residual distribution for m 2 <C mf 2 for the CfCa terms is very small compared to the 
size of the non-global terms. We expect that power corrections in m^/Q 2 account for the 
remaining tiny discrepancy from zero. 

The nice agreement with EVENT2 for the CfCa color structure is a very strong check 
on our calculation. The m 2 dependence of the CfCa distribution has been fully accounted 
for by our calculation of the non-global terms, up to a very small correction likely arising 
from power corrections in m^/Q 2 . 
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Figure 7: The 0(a;?), CpTRnj distribution from EVENT2 with various sets of terms removed, 
as in Fig. § The upper left plot shows the distribution directly from EVENT2. The upper right 
plot has the known global terms subtracted. The lower left plot has the single log in the cumulant 
subtracted. The lower right plot has the global terms and all the non-global terms from our 
calculation subtracted. The red line marks the value of m\ used. Note the changing vertical scale 
between the plots. The high region of the distribution has large numerical errors; see the discussion 
in the text for more details. 



In Fig. [?] we show the distributions in Eq. ( |5.6| ) for the CfTruj color structure. Note 
that there is no double log term to subtract for this color structure, so the lower left plot 
subtracts the single log. Recall that the parameter choice NPOW1 = NPOW2 = 4, needed 
for reliability in the small m| region, has introduced large statistical and systematic errors in 
the large m?, region, to the point that there is a clear deviation of the EVENT2 results away 
from the true distribution for m^/Q 2 > 10 -3 . However, we still see consistency between 
our calculation and the EVENT2 results in the small m| region. These uncertainties at 
large m| make the result less definitive as the CfCa color structure, but it nonetheless 
gives us confidence in our results. 

5.2 Logarithms in the Hard Regime of One Hemisphere 

We have seen that our calculation of the soft non-global logarithms matches well to the 
results of EVENT2. These soft non- global logs appear in the region of phase space where 
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both jet masses are small, mf 2 ^ Q 2 j an d hence showed up as elements of the leading- 
order dijet factorization theorem. When one hemisphere jet mass is large, of order the hard 
scale, 

ml < mf ~ Q 2 , (5.7) 

the structure of the logarithmic terms in m-i is different. In this regime, the left hemisphere 
has two or more energetic jets that give a large value of m\ and the right hemisphere has 
one collimated jet that gives a small value of 777,2. From the point of view of the dijet 



factorization theorem in Eq. (1.4) the new logarithmic terms that contribute for the region 
in Eq. fl5.7p are power-suppressed terms, of the form (m 2 /Q 2 y ^ln fe (m|/Q 2 ) with j > l. 7 
We will refer to these terms as coming from the hard regime of m\. 

The logarithms of 777.2 in this hard regime have both global and non-global sources. 
Global logarithms come from collinear emissions from the jet in the right hemisphere and 
soft emissions from all energetic jets. Non- global logarithms come from secondary soft 
emissions from soft gluons, as in our calculation of non-global logarithms in the dijet 
regime. These non-global logarithms will be different than in the dijet regime, and require 
additional calculations. Since non-global logarithms in the hard regime start at 0{p? s ) we 
will only consider global logarithms from the hard regime here. 

To validate the above picture we will consider 0(a 2 ) logarithms in the hard regime 
and compare them to the results from EVENT2. Specifically, we will consider the double 
cumulant distribution with m\ ~ Q and 7772 <C Q- The soft non-global logarithms that we 
have calculated dominate in the small m\ region, and are only part of the contributions 
to the double cumulant when 777^ ~ Q. For this hard regime in 777^ we will compute the 
double logarithmic dependence on m c 2 for the CfCa color structure. 

At 0(a 2 ), the hard regime logarithms come from collinear and soft emissions from a 
configuration of three hard partons (a quark, anti-quark, and gluon). To have a large value 
of mi, two of these partons must be in the left hemisphere. The tree- level cross section for 
this 3-parton configuration is the well-known 



da 3 _ a s C F x 2 + x q 



2 1 „2 

dx q dXq 2TT (1 — X q )(l — Xq) ' 

where x qt q = 2E q> q/Q are the energy fractions of the quark and anti-quark, and x g = 
2E g /Q = 2 — x q — Xq is the energy fraction of the gluon. To determine the double logarithms 
in 7772, we consider an emission of a soft gluon from any of the 3 jets or a collinear gluon 
from the jet in the right hemisphere. For the double 7772 logarithm the sum of soft emissions 
depends only on the color of the parton in the right hemisphere. Thus the sum of all soft 
and collinear contributions yields a multiplicative factor to Eq. ( |5.8[ ), 

a, „ , /777 c2 



where Cj is the Casimir invariant of the jet in the right hemisphere {Cf for quarks or 



7 Using SCET in the dijet limit these terms are described by factorization theorems involving power- 
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Figure 8: The EVENT2 distribution for the CfCa color structure in m^/Q 2 before and after 
subtracting the double logarithm of m^/Q 2 from the regime of hard mf given in Eq. (5.10), for 
mf/Q 2 = 0.2. In both distributions we have subtracted the global dijet terms and the soft non- 
global terms. After we subtract the additional double logarithm coming from the regime of hard 



,<'2 



of: 



, the remaining distribution at small m<i is flat, and therefore contains only a single logarithm 

7Q 2 - 



anti-quarks, Ca for gluons). Therefore for the CfCa double logarithm the gluon jet is in 
the right hemisphere, and quark and antiquark jets are in the left hemisphere. All that 
remains is to convert the distribution differential in x q and x q in Eq. ( |5.cf ) to a cumulant 
in m\ . Up to power corrections in m 2 /Q 2 we can use m\ = Q 2 {x q + x q — 1) to determine 

x u 



The gluon being alone in a hemisphere requires 2 — x„ 



> Xq,q which yields 



the limits on x q . We also multiply by 1/2 since we are requiring that the left hemisphere 
mass be larger than the right, which restricts the gluon's angular phase space. Thus in the 
hard regime for m^ 2 , the CfCa double logarithm in m^ 2 is 



S hard (mf ,mf) 



c2 



,c2\ 



a a 2 s C F C A 1 r 1 
Q 4 (2vr) 2 2 



dm 2 



o 



cro ( a s \ 2 C F C A mf 
Q 2 



2tt 




2\2 



x 2 + (l-x q + m 2 /Q 
(l-x q ){x q -ml/Q 2 ) 



In 2 



m 



c2 



Q' 2 



,c2 



c2 



3 2>m\" / m 1 
2 + ~Q 2 ~ + V 2 + Q2 



ln 



»7 



c2 



+ 



Q 



2 r 



<-2 



m 



c2 



4Li 2 



/2m ( 



c2 



Q 2 

f 2Li 2 I 



+ 41n 1 



m 



c2 



Q 2 

(Amf{Q 2 - mf) 

'Q 1 



In 



<2 



Q 2 



Q 2 
+ 4Li 2 



-2mf 



in 



c2 



Q 2 



m 



c2 



Q 2 



(5.10) 



suppressed hard, jet, and soft functions, order-by-order in j (see e.g. [^2[). When N jets in the left 
hemisphere contribute to mi, these terms can also be described by an (N + l)-jet factorization theorem in 
SCET (see e.g. @ g g). 



- 36 - 



The overall m'f/Q 2 shows that this result is power suppressed for m^ 2 <C Q 2 , and in this 

limit the terms in {• • • } reduce to — 2 ln(m^ 2 /Q 2 ). If we set m^ 2 /Q 2 to the kinematic upper 

limit of 1/3, then Eq. ( |5.10|) reduces to the result in Eq. (|L~9| ) with m^/Q 2 = Pr. 

We can compare this additional double logarithm of m?? /Q 2 from the regime of hard 

mf 2 to EVENT2 in the same way that we compared the soft non-global logarithms. In Fig. || 

we show the effect of subtracting the CfCa double log from the EVENT2 distribution for 

the mixed distribution-cumulant cross section, da /dm 2 (ml 2 ), computed as in Eq. (^). 

Choosing a value of m^ 2 = 0.2 Q 2 , we first subtract the dijet global logarithms and soft 

non-global terms in m^ 2 given by Eq. ( |3.36| ), and plot the result as the red (upper) points 

in Fig. |8| We then subtract the double logarithm in m^ 2 from the hard mf 2 regime given in 

Eq. ( |5.10| ) and plot the result as the black (lower) points in the figure. It is clear that after 

subtracting the double logarithm the remaining distribution contains just single logarithms 

of mf, and therefore Eq. ( |5.10| ) correctly accounts for double logarithms from the hard 
m c2 re gi me _ 



6. Conclusions, Simple Generalizations, and Outlook 

We have calculated the dijet hemisphere soft function to order a 2 in both position and 
cumulant momentum space, S(x\,X2, p) and S c (£\,l\,\x). In doing so, we uncovered the 
full non-global structure of these functions, that is, the dependence on the ratios x\jx<i or 
which is not fixed by renormalization group invariance of the factorized dijet cross 



section Eq. (1.4). We re-derived the non-global double logarithm whose presence in related 



observables was pointed out by p0[ 21], obtained new results for the single logarithmic 
terms, as well as determining the full non-logarithmic non-global structure in the dijet soft 
function and related observables. Eqs. ( 3.3C| ) and ( 3.36| ) are the main results of this paper. 



These results go significantly beyond previous discussion of non-global structures in jet 
cross sections at 0(a 2 s ). 

The discovery of these non-global structures resolves the question of the structure of 



the dijet soft function. The ansatz of Ref. [42| that only double logarithms of x\jx2 and a 



constant could exist in S(x\, X2, p) due in part to x\ «-)■ X2 symmetry, is incomplete. Many 
other functions of X\jx2 appear, including a single logarithm ln(xi/x2 + ^2/^1) that is 
symmetric under x\ 0x2- 

A key piece of our strategy to simplify calculation of the dijet soft function was to 
take advantage of renormalization group properties and infrared finiteness of the full soft 
function S(x\, X2, p) or S(£±, £2, /u) so that all the non-global dependence at 0(a 2 ) could be 
extracted only from those diagrams in which two particles in the final state enter opposite 
hemispheres of the dijet event. This procedure determines all the terms in the soft function 
except for the constant term, which we take from a recent calculation of this constant in 



the 0{a 2 s ) thrust distribution by |4C]. Cancellation of IR divergences between same and op- 



posite hemisphere final-state configurations leaves over logs of x\jx2 or £i/£2, illuminating 



how NGLs find their origin in the structure of a soft function [48|. 

In this paper we explored non-global logarithms that arise due to the sensitivity of 
an observable to soft radiation at two different scales in sharply divided regions, i.e. m\ 
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and m,2 in opposite hemispheres. These observables can naturally be taken to be similar, 
mi ~ m 2, or parametrically disparate, m\ 3> ni2. The NGLs of ln(mi/m2) become large 
in the latter case, and are distinguished from global logs ln(mj/Q) which are the only large 
logs in the former case. Ref. [^] distinguished observables with soft radiation measured 
everywhere, from those in which soft radiation is probed with a single scale in only a 
portion of phase space, e.g. m\ in one hemisphere (the NGLs originally studied in pCp ). 
We found in the case of hemisphere masses that NGLs of the latter type can be derived 
from those where soft radiation is measured everywhere, by integrating the variables for 
the unmeasured regions up to the hard regime. Thus from an effective field theory point 
of view which focuses on hierarchies of scales, we find it more economical to think of NGLs 
as being of a single type. We consider the measurement of independent observables like 
m\ and ni2 that cover all of phase space, but separately covering sharply divided phase 
space regions. These kinematic variables can naturally be taken similar (mi ~ 1TI2), or 
can be taken as disparate (mi 3> 7712). In the former case it suffices to treat NGLs in 
fixed order perturbation theory. In the latter case we have large NGLs that appear from 
soft radiation that probes the sharply divided regions, and NGL resummation [20] must be 
carried out. Effectively we consider NGLs with all variables unintegrated and derive NGLs 
for other cases by integration. Keeping the variables unintegrated facilitates the task of 
distinguishing logs of global and non-global origin. 

We performed several cross-checks of our results, by considering general covariant 
gauge, by projecting the dijet soft function onto various observables and comparing to 
EVENT2 predictions for those observables. These comparisons revealed existing EVENT2 
extractions for the difference between heavy jet mass and thrust soft function constants to 
be in excellent agreement with our analytic result. We also compared our prediction for 
the double hemisphere mass distribution to our own EVENT2 runs, and again found good 
agreement. Finally we considered the regime of one hemisphere mass mi becoming hard and 
generating additional global logs of 777-2 /Q- Our calculation of the additional double log in 
this regime also matched our EVENT2 runs very well. We also verified the above-mentioned 
relation between observables containing NGLs with soft radition measured everywhere, 
versus measured only in a region of phase space, by integrating 7771 all the way up to its 
kinematic limit and reproducing the global and non- global double logs of pr in Eqs. 
and Q at 0(a 2 s ). 

State-of-the-art fits for a s from e + e~ — > jet data currently rely on N 3 LL resummations 
of event shape distributions fl6|, |l^, and thus on the fixed-order 0(a? s ) dijet soft function. 
Our new results make possible improved extractions of the strong coupling by removing a 
source of uncertainty in the perturbative calculation. 

The 0(a 2 s ) soft function in position space, in particular, makes possible a fully analytic 
N 3 LL resummation of the doubly differential cross section Eq. (1.11) in position space 
when X1/X2 ~ 1, and is also a crucial ingredient in the N 3 LL resummation of the doubly 



differential momentum space distribution Eq. ( |l.4|) or cumulant Eq. (1.5). 

One simple generalization of our results is to the case of hadron colliders. For pp 
and pp collisions an incoming dijet hemisphere soft function appears for the event shape 
called beam thrust tb (or 0-jettiness), which in the dijet limit tb <C 1 provides a veto on 
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central jets. The only difference between our soft function in Eq. ( |2.6[) and the incoming 
one are that the Wilson lines are incoming from — oo rather than outgoing to +00. This 
reflects itself in the signs of i0 + factors in the eikonal denominators in the momentum space 
amplitudes. Since for our opposite hemisphere calculation the signs of those i0 + factors 
drop out, our analytic results for the x\/x2 and if/i^ dependence immediately carries over 



to the incoming case. In Ref. [55] it was proven that the //-dependent logarithms are also 
the same for the outgoing and incoming dijet hemisphere soft functions. Hence at 0(a 2 s ) 
the only thing that can differ between the two cases is the thrust constant sr? wGA,nt \ (\i 
would be interesting to compute these constants directly for the incoming case.) Thus our 
results already make possible an improvement of predictions for jet observables in hadron 
collisions. 

Another simple generalization of our results are to Wilson lines in different color repre- 
sentations, such as octets for incoming or outgoing gluon dijets. It is straightforward to see 
that our computations obey Casimir scaling, so we can simply modify the CVs to get re- 
sults for other representations. For the Q-,%-, and Q graphs only traces that obey Casimir 
scaling appear. For the I and T graphs non-abelian exponentiation and the generator 
commutation relations guarantee that Casimir scaling also applies. 

Looking ahead, relating non- global logarithms to the soft function in a factorization 
theorem in the context of effective field theory as we have done also opens the door to the 
development of renormalization group techniques to resum non-global logarithms. This 
would make it possible to go beyond the resummation of NGLs in the conjectured form 
given in |2C] which is parameterized in terms of coefficients determined numerically in the 



large- Nq limit. It is important to keep in mind that such NGLs (and other non- global 
structures as we found here) will appear not only in the single or double-hemisphere mass 
distributions, but also in most jet cross sections in which jets are defined with one of the 
well known jet algorithms, which leave the regions inside and outside these jets sensitive 



to different soft scales [31 



Note added: While this paper was being finalized, Ref. |)(| appeared calculating the 
dijet hemisphere soft function in momentum space, S c {l\,t\,ii). Our result in Eq. ( |3.36j ) 
exactly agrees numerically with theirs, even though the two have different algebraic forms. 
In our algebraic result the logarithmic singularities are separated out and the remaining 
polylog dependent contributions [our /(^f/^f ) + f {^-2/^-1 ) given by Eq. ( |3.37| )] are bounded. 
An analog of our two dimensional position space result S(xi, 22 > /") in Eq. ( |3.30 ) was not 
obtained in Ref. |5(|. On the other hand, unlike our paper Ref. |5(| does in addition 



compute analytically the thrust constants, finding = — — fg|vr 2 + ^vr 4 + 

and = |t + 2T 71 " 2 — T'' 3 ' These agree numerically with the first four digits of the 
result from Ref. |4^] quoted in Eq. (|2.22| ). In the publication corresponding to [pCj] , which 
is Ref. 57], analytic results were also obtained for the thrust constants, and these agree 
with those of Ref. [56]. 
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A. Anomalous Dimensions 



In Sec. 2.2 we reviewed how the //-dependent terms in the soft function can be constrained 
using renormalization group invariance of the cross section Eq. (1.4). Here we record the 
soft anomalous dimensions that are needed to carry this out explicitly to G(o? s ). 



The exponents K, uj appearing in the evolution factors Ug(x, (j,, [Aq) in Eq. ( p. 13 ) are 
defined by 

f 1 * du' ( , fi' , \ 

K(T cusp ,j s ,^^o) = —r -2r cusp [a: s (//)]m h 7s[a« (//)] , (A. la) 

J no M V Mo ) 

w(r cusp ,/i,/i ) = -2 / — r cusp [a s (//)] , (A.lb) 
J no A 4 

where the cusp and non-cusp parts of the anomalous dimension r cusp and 75 are defined 
by Eq. ( p. 11 ). We define the expansions of the anomalous dimensions in a s as 



00 k+i 00 fe+i 

rw« s ] = £© r c fc usp , 7s[« s ] = £(^) tS- (a.2) 



fc=0 fc=0 



At Q(a s ) and 0(aj) M, 59], the cusp anomalous dimension is given by 



C, = 4C F (A.3a) 

■«..!^-C^, (A.3b) 

and the non-cusp anomalous dimension is given by 

l°s = (A.4a) 

4 = C F C A (_f + If! + 28C3 ) + C F T R n f (f - ^) . (A.4b) 
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We also need the coefficients of the beta function, 



(3[a s 



k=0 



Utt/ 



(A.5) 



where 



A) = ^-C A - ^T R n f , 



3 

34 

y 



3 

—C a Trti( - AC F T R n F . 
3 



(A.6a) 
(A.6b) 

In terms of these beta function coefficients, the two-loop running coupling is given by 
1 



! , A) / m A*' 
1 + — a s (/i)ln — 

Z7T // 



(A.7) 



These are all the quantities needed to determine the /^-dependent pieces of the soft function 
to 0(a 2 s ). 8 

The quantities K, uj to 0(a 2 s ) written out explicitly are 



7s,M,Mo) 



«s(m) 



47T 



Mo Mo 



(A. 



+ 



w(r cusp , fi, no) 



Air 
«s(m) 

4"7T 



jr° usp /3o In 3 £ + ( 7 §ft - T^p) In 2 ^ + 7 | In ^ 



Mo 



Mo 



-2r° In — 

z,± cusp 111 , 
MO 



Mo. 
(A.8b) 



These expressions can be used to determine the parts of the soft function dependent on 
logs of |U, Eq. ( 3.35D in position space and Eq. ( |3.43 ) for the double cumulant in momentum 
space. 



B. 0(a 2 s ) Diagram Results for the Opposite Hemisphere Soft Function 

Here we present the details of our calculation of the 0(a 2 ) opposite hemisphere soft func- 
tion, in both position and momentum space. The opposite hemisphere soft function in 
momentum space is given by 

J^J^Mh,k 2 )M^%(h,£2)C(k 1 )C(k 2 ), (B.l) 



To achieve NNLL accuracy in the resummed soft function Eq. (2.12) it is actually necessary to know 
the cusp anomalous dimension and beta function at 0(0:3) as well, but in this paper we are only concerned 
with the fixed-order result. 
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where C(k) is the cut propagator, Eq. (|3.5| ), we sum over the squared matrix elements from 
different classes of diagrams Aj where j = {I, T,Q,7i, Q} (see Fig. [I]), and the momentum 
space measurement function is 

M l £%(£i,t2) = - K)8(£ 2 - k+)9(kf - k^)9{k^ - k$) 

+ 5(ti - k^)5(£ 2 - k+)6(k^ - k+)9(k+ - fcj ) . 

Below we will denote the integrated contribution of each diagram j to S opp in Eq. ( [B.lD by 
the symbols I, T, G,T~l, Q themselves. The position space opposite hemisphere soft func- 
tion is given by replacing the measurement function above by its position space analog, 
A4\^ 2 (xi,x 2 ) in Eq. fl3.12| ). The integrated result analogous to Eq. (BA) for each di- 



agram's contribution to S opp (xi, x 2 ) will be denoted I,T,G Q- More generally, the 
measurement function for other observables, Mkxkn can ^ e use d with our results for the 
double cut squared matrix elements, Aj, to obtain the contribution to the soft function for 
other observables. 

The observable and hemisphere geometry determine some natural variables for the two- 
loop calculation. The measurement function and the diagrams are symmetric in k\ ■<-)■ k 2 , 
and therefore we focus only on the case where k\ contributes to i\ and k 2 contributes to 
t 2 . The initial matrix elements are straightforward to write in light-cone coordinates, and 
the measurement function A^j^^ in Eq. fl3.12| ) implements the conversion 

k^-tti, k+ -> l 2 . (B.2) 

It is natural to then scale out by these variables. We find that a pair of variables particularly 
useful for the calculations are 




/^l / k-^ / k-^ / k-^ j k-^ £ 



u=J^2_ /^I?.. (B.3) 



k-^ k<2 y ^2 y k^ J k<2 y k<£ 



The phase space cuts constrain z > 1, u > 0. Physically, the ratios kf /k± and k^ /k 2 are 
related to the polar angles of the partons with respect to their jet directions. For example, 
kf /ki = tan 2 (6*i /2), where 9\ is the polar angle of the parton with respect to the jet 
(thrust) axis. 

The z variable measures how close to the boundary the combination of the two partons 
are: as z — >■ 1, both partons move towards the hemisphere boundary, and as z — > 00, the 
angular separation between the partons grows, u measures the relative angles of the two 
partons. It is worth noting that divergences from z — > 1 play a role in the appearance of 
non-global logarithms. 

With the exception of the I matrix elements that are straightforward to evaluate, we 
find a common set of transformations are useful to simplify the integral form of the matrix 
elements. The matrix elements depend on the set of variables 

{kf,ki,k+,kz,<j>}, (B.4) 

where <j) is the angle between k\ and k 2 in the plane transverse to the thrust axis. Our 
matrix elements depend only on cos 4>, which arises from k^-k^ ■ Starting with the two- loop 
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measure, we can perform some of the integrals to obtain 
d d h d d k 2 ... ... 1 (e^) 2e 



/d k\ d k 2 1 (e 7E ) f°° 

(2ir) d (27rp C ^^ C ^ 2 ) = (I6vr 2 ) 2 T(l - e) 2 J ^ 1 ^ 2 ^ 2 ^ ^ ^ 2 ^ 2 



V/ 2 r(i 



r(i - e) 



-1 

/*7T 

/ ^sin" 2e <}>. (B.5) 
Jo 



The measurement function then allows us to perform the k l and k 2 integrals. For the 
opposite hemisphere contributions shifting to z and u requires the change of variables 

dkf dk 2 + ,_ WI _ ^ ,+s _ f°° , f z , 2z 



f°° f z 2z 
6(k+ > k^)9{k^ >k+)= dz du— . (B.6) 

Jl Jl/z U 



o ^1 k 2 Ji Ji/ Z 

This change of variables is convenient because in all the diagrams the u integration is 
independent of e and can be easily performed. Additionally, the 4> dependent terms in the 
matrix elements, which give hyper geometric functions when integrated, are independent of 
u. This means that we can simultaneously perform the u and </> integrals, leaving only a 
single integral over z to be performed in each diagram. Standard techniques can then be 
used to evaluate these z integrals. 

In addition to a common strategy for evaluating the integrals, several common func- 
tions arise in the calculation. Using the z, u variables, a function of r that appears in nearly 
all the diagrams is 

g(z, r) = In ( . ( 1 + r / z ) . ( B .7) 
V ; \{z + r){\ + zr) ) v ; 

This function is symmetric in r — > 1/r and scales as r for small r. In addition, there are 
two common hyper geometric functions that arise from integrations over </>: 

(^l 2 T(\ -e)\ r g . n _ 2e 1 = _1 

J r Y 1 + z 2 -2zcos<j) z 2 + V v ' K ' 



7T 



r(i - e) 

'It) I # 5in ' 2 ' (i + g ^ = -<?w ^- 2 > - 2A '- 2 )i • 



where 



'1 4z 
fi(z) = 2 F 1 ( — , 1, 1 — e, (1 + z)2 

/ 2 (z) = ^(1,2,1-6,^-^) . (B.9) 



It is necessary to expand these functions to 0(e 2 ): 

h(z) = {l + z) z 2e (z - l)- 1 -^ [1 + 2e 2 Li 2 (l/z)] + 0(e 3 ) , 



2e iA-3-2e 



Mz) = (l + z)z" (z-1) 



l + z 2 + 2e(l + z) 



- 2e V - 1) In [—\ + 2e 2 (l + z 2 ) Li 2 / ' 



+ C(e 3 ) (B.10) 
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Note that we have summed up certain terms at higher orders in e in these expansions. 
Throughout the calculation we made use of the HypExp package [6C, 31 1. The divergent 
structure in e of the vacuum polarization diagrams is quite complex, and careful treatment is 
necessary to ensure the divergences are regulated and terms are not forgotten. The difficulty 
in evaluating these diagrams comes from the hypergeometrics; using certain hypergeometric 
forms we can see more clearly the regulation of the divergences in the integrand about z = 1 . 

B.l Form of the Matrix Elements in Momentum and Position Space 

Before giving results for the matrix elements, we discuss the general form of the soft 
function. We determine the matrix elements in both momentum and position space. In 
position space the soft function is a simple function of the variables x\ and X2, which are 
conjugate to £\,£2- In momentum space the £\,£2 dependence gives distributions (rather 
than functions), but the result does not take its simplest form in these variables. Instead, 
we use variables 



£ 2 f°° ,. f°° , , 1 



, ' IT 



dsdr — . (B.ll) 



The distributions in terms of these variables are simple: all matrix elements contain a dis- 
tribution of s and non-singular functions of r. The general form of the opposite hemisphere 
terms in momentum space is 

S opp (£ 1 ,£2) = ^S opp M(£i,£ 2 ), (B.12) 

i 

where 

S opp ®(£ 1 ,£ 2 ) = ACiS- 1 ' 26 ^) (B.13) 

and Ci is a color structure, (Cp, CfCa, or CpT^nf), F(r) depends on the class of diagram, 
(X,T,g,H,Q), and in MS 



/O.N2 ( e 7y )2e 

~ \2ir) r(l-e) 2 ' 
The transformation from s,r space to the position-space matrix element is 



(B.14) 



5 opp W(xi,x 2 ) = Ad / dsdr — e-^^/^+^h- 1 -^ Fir). (B.15) 
Jo 2r 

Convergence of the Fourier transform requires x\ 2 to have a small negative imaginary 
component: 

xi-iO + , x 2 -i0 + . (B.16) 

Because both s and r are positive, the sign of i0 + relative to x\^ is unchanged for the 
entire Fourier transform. We will often write just x\ 2 with the prescription in Eq. ( |B.1G| ) 
implied. 

In all matrix elements the function J-(r) is non-singular, meaning it can be split into 
a pure constant plus a nontrivial function of r, 

T(r) = F + F 1 (r) , (B.17) 
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where F\(r) is symmetric in r — > 1/r, and vanishes as r — > 0,oo. The s integral can be 
done immediately, giving 



S opp %i,x 2 ) =AC i T(-ie) 



dr 
o r 



[Fo + F^r)] 



Ac 



The r integral in the Fq piece can also be done immediately, using 



00 dr 
r 



Xl 



+ l2V r 



•Ic 



\IX1IX2) 



\2c 



r(-2e) 2 



Then each contribution to the opposite hemisphere soft function can be written 



S opp M(xi,x 2 ) = ACi(i Xl ix 2 ) 2e < 



00 dr 

r(-2e)^F + r(-4e) / — F x (r) 
r 



b jr 
r + Vb 



(B.18) 



(B.19) 



ic 



where 



X\ — 1 



0+ 



(B.20) 

(B.21) 

x 2 - i0+ 

The function i*i(V) always vanishes sufficiently quickly as r — > and r — > 00 so that the r 
integral introduces no additional powers of 1/e, and thus the integrand can be expanded 
and truncated at 0(e). Splitting the integral into the regions < r < 1 and 1 < r < 00, 
and using the symmetry of F\{r) in r — > 1/r, we can rewrite the last integral as 



f 1 dr„ , \ 




/ -Mr) 


2 + 4eln f r + - 


10 r 





(B.22) 



The last integral displays logarithmic dependence on b and 1/6, which can be extracted by 



ln(r + - + 6+ i 

r 



In ( 2 + 6+ - J + In 



. 2+6+i 



(B.23) 



Then Eq. ( ^20| ) can be written 



S opp W(xi,x 2 ) = AC i (ix 1 ix 2 ) 2t < r(-2e)^F + r(-4e) 



2 + 4e In 2 + b + 



+ 4eT(-4e) / — Fi(r)bi 



. 2+6+i 



(B.24) 



In each term, the i*b,i also have expansions in e, and should be kept/truncated to the 
appropriate order. In particular, we must keep terms of up to order e 2 in Fq, but only 
to order e in F\(r) (and only e° in the last integral). Note that in the last integral the 
integrand now vanishes at both endpoints, r — > (due to F\(r)) and r — > 1 (due to the 
log). 

We now give the result for each group of diagrams in both position and momentum 
space, as well as the intermediate step where only the z integral remains. 



-45 - 



B.2 X and T diagrams 

The matrix element for the independent emission diagrams in Fig. |l](a) is a sum of contri- 
butions with color factors C F and C f Ca, Ax = A x c 2 + Ax,c F c A - The Cj, piece is 



l Q„4.,4e/-i2 



8<?V e C 2 



F 

1 1 

A/-j^ »VQ 



(B.25) 



This diagram is the product of one-loop diagrams with one gluon in each hemisphere 
(neglecting the delta function in the one- loop diagrams for the hemisphere with no parton) . 
In position space, the result of the diagram's contribution to S opp is 

X C 2 = ylCj,(^iix 2 ) 2e r(-2e) 2 -i . 
The result in momentum space is 

X c% = ACfihh)- 1 - 2 ^ . 
The C f Ca contribution to the independent emission matrix element is 

4e (^l + ^2X^1 + k 2 )+kfk 1 + k 2 k 2 



(B.26) 



(B.27) 



Ax,c F c A = ~ 2 9 CfCaU 



{k\ + k 2 ){k 1 + k 2 )kfk 1 k 2 k 2 



(B.28) 



To evaluate this diagram we use the variables u = kf jk 2 and v = k 2 /k l . An intermediate 
step in the calculation is 



%C F C A = —ACfCas 



-l-2e 



1 



+ 



l/r 



du \ dv (uv)' 



U + V 



(1 + «)(! + «) 



(B.29) 



These integrals can be evaluated exactly in e. The leading divergences in Ic F c A cance l 
with the T diagrams, so we present the results for the sum of the diagrams. 

The matrix element for the single 3-gluon vertex T diagrams in Fig. |l|(b) and (c) is 



A T = 2g 4 C F C A ^ € 



h~^~ h~ -X- h~ 
1 2 1 2 

2h ■ k 2 



1 + 



1 1 2 2 



1 



U+ j.— U+ u- 



(B.30) 



(kt + k^)(k^ + k 2 )_ 

Using the variables and transformations described in Eqs. ( p3.2[) - (|B.10|) , the intermediate 
step with the remaining z integral in the calculation for the diagrams' contribution to S opp 
is 

z + 1 



T = AC F C A s 



-1-2<E 



dz z 



-1-e 



21nz + 



z- 1 



g(z,r) 



(B.31) 



Since the leading divergences in the Xq f c a an d T diagrams cancel, we sum them. In 
position space, the result for their contribution to S opp is 



+ 4r(-4e) 



7T 



Zc F c A +T = AC F C A {%xx%x^\ r(-2e) 2 ( — + 4( 3 e + ^ 



2^ 



14vr 4 



7T 4 vr 2 M 



+ 



2tt 2 



r 3 2e C3 H 

X 2 Xi ' 



X 2 J \Xl 



(B.32) 
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where the last two terms are given by the integral 
■ i 



41n(l + r) 41n(l + i) , / N / 1 
; + ^ V , rJ - 2 ln(l + r) In 1 + - 



1 + r 



1 + 



7rr(b) + J=rr(l/b) 



r 

(B.33) 



where 



7T 



In 4 b 



2^ 



7T 



^rr(&) = ~ + + ( 2Ca - =jM In(l + 6) - ln 2 Mn(l - 6) - y Li 2 (l - 6) 



41n6Li 2 (&) + (Li 2 (l-6)) 2 + 6Li 3 (5) + 21n6Li 3 (l - b) . 



(B.34) 



The sum of the two terms Fixib) = J~xT$) + J~XT(^/b) is a finite function of b, in fact 
vanishing as b — > or b — > oo. 

In momentum space, the result for the Ic F C A ,T diagrams' contribution to S opp can 
be written 



Z Cf c a +T = AC f C a s 



-l-2e 



+ 4C 3 e + + £rr(r) + £rrO-/r) 



3 ' 45 
AC7 J? C7 A s- 1 - 2e [i ; b + i i i(r)] 



(B.35) 



where 



7"zr(r) = ln(l + r) In 1 + 



1 



1 + r 



ln(l + r) 



+ I - In 2 r 



lnr + 



4r 
1 + r 



1 + r 
1 - r 

! lnr - 4 ) ln(l + r) 

1 + r 6 



+ 2 lnr Li 2 (-r) 



7T 



(B.36) 



The r-dependent functions appearing in Eq. (B.35) are the Ic F C A + T contribution to the 



function F\(r) defined in Eq. ( p3. 17| ) . This contribution is explicitly symmetric in r — > ~ 
and vanishes as r — > or r — > oo, as advertised. This is also true of the other diagrams, 
which have the same form. 



B.3 Q and % diagrams 

The matrix element for the gluon loop Q diagrams in Fig. |](d) is 

1 



Ag = g C F C A V 



k 



(2fci • k 2 ) 2 (k+ + k+)*(ki +k 2 y 

+ kt)(K +k 2 )[- I6h ■ k 2 - 2(1 - e)(kf - k+)(k^ - k 2 ) + 4(A:+ + k+)(k^ + 
+ k+k 2 + fcjf + (fcr + A^) 2 [(1 - e)(fc+ - k+) 2 - 2(k+ + k+) 2 - k+k+} 
+ (fe+ + k+f [(1 - e)(A:r - k 2 f - 2fc + k 2 ) 2 - k^k 2 ] } . (B.37) 
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The matrix element for the ghost loop H diagrams in Fig. |l](e) is 

1 



k 



An - g C F C A » hm + kt)2{k _ + k _ )2 

X *^ ^2 ) (^1 — ^~ ^2 ) (^1 ^2 — ^" ^1 ^2 ) — ^~ ^1 ^2 — ^" ^2 ) — ^" ^1 ^2 (^1 — ^~ ^2 ) } " 

(B.38) 

The ghost loop diagram cancels a set of terms in the gluon loop diagram, and so we evaluate 
them together. Using the variables and transformations described in Eqs. ( |B.2| ) - (jB.lCj ), 
the intermediate step with the remaining z integral in the calculation of the diagrams' 
contribution to S opp is 



G + H = A CfCas- 1 -^ dz z-*L—L^ g(z, r)h{z) 



+ 2(l-e)^—g(z,r)[f 1 (z)-2f 2 (z)] 
z z — 1 

+ - ^ (r + ,; i + ra) IAW - *AM) i ■ (B-39) 
In position space, the result for the diagrams' contribution to S opp is 



3 9 y ' V 9 54 ' 3' 3 

^M 2+ S + S) + MSH«(t)}' < B -» 

where the last two terms are given by the integral 

Jo r \ 2 + 6+ i j 3(l+rf 

2r(l + r) + (5 + 9r + 6r 2 ) ln(l + r) + (6r + 9r 2 + 5r 3 ) In f 1 + - 
= J r 6w(&)+^(l/6) J 
with 

- ^ln6Li 2 (6) +5Li 3 (6). (B.42) 
In momentum space, the diagrams' contribution to S opp is 

+ H = AC F C A s-^^Pguir) + Pguil/r)) (B.43) 
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where 



3(1 + r) 3 
2f 



3(1 + r) 3 



2r + (5 + 9r + 6r 2 ) ln(l + r) 



16 5vT 2 x o /31 9 \ , 

yr + — (1 + r) 3 - f y + 23r + 8r 2 j ln(l + r) 



(B.44) 



1 



(5 + 3r - 3r 2 - 5r 3 ) ln 2 (l + r) + (6r + 9r 2 + 5r 3 ) Li 2 (-r) 



In this case Eq. ( B.43| ) takes the form given by Eqs. ( EDI ) and J51?D with F = and 
B.4 Q diagrams 

The matrix element for the Q diagrams in Fig. |l](f) is 

1 



Aq 



C F T R nf fi 



Ac 



(2fci • A: 2 ) 2 (A;+ + A;+) 2 (fc 1 - + A^) 2 

X ^2 — ^~ ^2 ) — ^" ^1 ^2 (^1 — ^2 ) 

+ (2fci • £; 2 - A^fc 2 - fc7"fc^) (ftf + k+)(k^ + fe^)] . (B.45) 

Using the variables and transformations described in Eqs. ( |B.2j ) - ( |B.10| ), an intermediate 
step in the calculation for the Q diagrams' contribution to 5' opp is 

1 



Q = AC F T R n fS - 1 - e dzz~'l 



> , g(z,r)[h(z)-f 2 (z)] 

z z — 1 

1 ' \h{z)-2f 2 (z)) 



z + 1 (r + z)(l + rz 

In position space, the result for the diagrams' contribution to S opp is 

4 8tt 2 



(B.46) 



Q = AC F T R n f {i Xl ix 2 ) 2e \ T(-4e) ( -- + () 



/4vr 2 



X\ x 2 



--- In 2 + - + - +T Q — + Tq 



\ 9 "! / \ .!■■> .rj 
where the last two terms are given by the integral 



17 16vr 2 4, . 

IT " — + 3 C<3) 



,r 2 



In 



' r + 1 + 6 + | 
2 + 6 + ± 



3(1 +r) 



(B.47) 



(B.48) 



2r(l + r) + (1 + r) 3 In 2 + r + - + (1 - r 3 ) In r 



?Q(P)+F Q {l/b) 



with 



F Q (6) ee 



3(6-1) 



In 6 



3(6 -iy 



ln 2 6 



2 _ 4tt 2 \ 

3 ~ ~9~ ) 



ln(l + 6) 



(B.49) 



+ \ In 2 6 ln(l - 6) + \ In 6 Li 2 (6) - 4 Li 3 (6) . 
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In momentum space, the diagrams' contribution to S opp is 

Q =AC F T R n fS - 1 -^[P Q {r) + f Q (l/r)] , 

where 

4 



+ 



3(1 + r) 3 
8e 

3(1 + r) 3 



[-2r + (1 + r) 3 ln(l + r) + lnr] 

-r + — (1 + r) 3 - -(1 + r) 3 ln 2 r 
3 24 V ; 8 V ; 



(B.50) 



(B.51) 



- + 7r + r 2 j ln(l + r) + ^(r 3 - 1) ( ln 2 (l + r) + Li 2 (-r)) 



Again the form of Eq. ( |B.50| ) is that given by Eqs. (|B.13[ ) and ( |B.17| ) with Fo = and 
Fi{r) = P Q {r)+P Q (l/r). 

B.5 Total Opposite Hemisphere Soft Functions 

Adding up the results in Eqs. ( p.32j ), (B.40), and ( B.47| ) the total opposite hemisphere soft 
function in position space is 



2c 



S opp (xi,x 2 ) = 

a s y (e^V) 
2 J T(l-e) 2 

+ 



(i Xl ix 2 ) 2t { T(-2ef 



C 2 F - + C F C A [^ r +4( 3 e + 



27T 2 

3 



147T 2 



2r(-4e) - In ( 2 + ^ + ^ 
x 2 X\ 

,, „ i 1 n ^ 2 ^ 20 67vr 2 4^ 4 22 



C F T R n f 



X-2 



+ Ci?CA Fc F C A — + FcfCa — 
V \xij \x 2 

+ FcpTRUf 



2 4^ 2 /34 32vr 2 8 

3 + ~~9~ + \ ~9 ~ 27 + 3^ 3 



+ 



C F T R n f (t Cf t r u s ( ^ 



(B.52) 



where 



J~C f Ca 



\xi 



which are given in Eqs. ( B.34 ) and ( B.42 ) and 



Fc'pTnnf 



( -L 



(B.53) 



(B.54) 



which is given in Eq. ( B.49 ). In the final result for the soft function S(x\,x 2 , fJ>) in position 
space given by Eqs. ( 3.30 ) and ( 3.31| ), we chose to replace the single log ln(2+xi/a;2+X2/xi) 
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in Eq. ( |B.52| ) with \n{x\/x2 + x%/x\) and shift the difference into the sum of J- functions, 



resulting in the -Fq,jv functions given by Eq. (3.32). This is so that the remaining terms 
Fq,n{xi/x2) + Fq,n{x2/x\) in Eq. ( 3.31| ) are bounded along the real b = X1/X2 axis, as 
illustrated in Fig. ||. Otherwise they would have a pole at b = —1. 

In momentum space adding up the contributions from Eqs. ( |B.35| ), ( |B.43| ), and ( |B.5C| ) 
the total opposite hemisphere soft function is 

S ™(4,4) = (^) 2 ^Jc|( W2 )— 1 



2n) r(l -e) 2 Vl 2 > 



+ CfCa 



Ar(r) + Jrr(l/r) + ^(r) + ^(l/r) 



+ C F T R n f 



?Q(r) + T Q (l/r) 



(B.55) 



where Fzt(t)i ^ r gu( r )i an d -^b( r ) are defined in Eqs. (B.36), (B.44), and QB.51J ), respec- 
tively. 
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